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ABSTRACT 


The  basic  compilation  and  analysis  of  the  seismic  data  during  the  International 
Seismic  Month  have  been  completed.  We  have  become  involved  in  the  develop¬ 
ment  of  management  systems  for  global  seismic  data  and  network  documenta¬ 
tion.  We  are  beginning  to  focus  our  attention  on  earthquake  source  mecha¬ 
nisms  and  the  determination  of  focal  depth;  some  initial  results  are  described. 
We  continue  to  evaluate  the  effects  of  earth  heterogeneity  both  on  surface  waves 
and  body  waves.  Development  oi  a  seismic  computer  terminal  system  for  the 
ARPA  network  is  progressing  satisfa"torily. 
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SUMMARY 


This  is  the  twenty-first  Semiannual  Technical  Summary  report  describing  the  activities  of 
Lincoln  Laboratory,  M.I.T.,  in  the  field  of  seismic  discrimination.  These  activities  involve  re¬ 
search  into  the  fundamental  seismological  problems  associated  with  the  detection,  location,  and 
identification  of  earthquakes  and  nuclear  explosions.  We  are  also  concerned  with  the  develop¬ 
ment  of  methoc.s  for  the  handling  and  analysis  of  large  quantities  of  global  seismic  data,  both 
from* the  pcir.t  of  view  of  data  management  system  design,  and  also  to  facilitate  the  optimum  ex¬ 
traction  of  scientific  information  from  high-quality  digital  data. 

A  major  achievement  during  the  current  research  period  has  been  the  completion  of  the  basic 
compilation  of  seismic  data  for  the  International  Seismic  Month  (ISM).  These  results  have  been 
published  in  the  form  of  two  Lincoln  Laboratory  Technical  Notes  (1974-14  and  1974-15),  and  the 
highlights  of  these  final  reports  are  described  in  this  SATS.  The  event  list  will  form  an  ideal 
data  base  for  a  variety  of  studies  in  seismic  discrimination,  and  the  experience  gained  during 
the  study  is  expects  i  to  have  a  significant  impact  on  any  future  operational  monitoring  scheme. 

We  are  becoming  heavily  involved  in  the  design  of  the  data  management  system  that  will 
accompany  the  new  seismic  instrumentation  currently  being  deployed  by  DARPA  NMRO.  Our 
basic  charge  is  to  ensure  that  the  acquisition,  flow,  and  storage  of  these  new  digital  data  meet 
the  requirements  of  those  groups  actively  engaged  in  studies  related  to  se.smic  discrimination. 

A  related  aspect  of  this  work  involves  the  exploitation  of  the  full  power  of  the  computational  facil¬ 
ities  of  the  ARPANET  for  the  storage,  transmittal,  and  retrieval  of  both  seismic  data  and  the 
associated  network  documentation.  Development  of  toese  systems  is  still  in  an  early  stage,  and 
a  brief  progress  report  is  included  in  this  SATS. 

With  the  completion  of  the  bulk  of  the  data  analysis  a.-*uociat<-d  with  the  ISM,  we  are  embark¬ 
ing  on  a  series  of  investigations  into  the  mechanism  and  focal  depth  of  the  seismic  source.  Some 
initial  results  are  outlined  in  this  SATS.  It  is  shown  that  average  Mg-mb  and  moment-Mg  re¬ 
lationships  place  bounds  on  acceptable  theories  for  the  spectral  characteristics  of  the  seismic 
source.  Information  about  these  relationships  can  be  deduced  from  frequency-magnitude  curves. 
Other  studies  are  concerned  with  the  effects  of  uppcr-n  antle  structure  on  focal  mechanism  solu¬ 
tions,  and  the  determination  cf  source  dimensions  for  earthquakes  in  a  localized  geographical 
region. 

Investigations  into  the  determination  of  the  focal  depth  of  seismic  events  are  proceeding  in 
several  directions,  two  of  which  are  described  here.  It  is  demonstrated  that  maximum  entropy 
spectral  analysis  may  be  applied  to  the  identification  of  depth  phases  by  the  location  of  cepstral 
peaks.  Tins  method  appears  promising.  Another  study  continues  our  analysis  of  the  estimation 
of  focal  depth  from  surface  wave  spectra.  It  is  shown  that  certain  earth  structures  that  are 
consistent  with  observed  surfacp  wave  dispersion  curves  can  lead  to  very  inaccurate  depth  esti¬ 
mates.  It  is  not  yet  clear  that  this  method  has  any  useful  application  to  the  seismic  discrimina¬ 
tion  problem,  particularly  at  low  magnitudes. 

Our  studies  of  the  effects  of  earth  structure  and  path  variations  have  continued,  and  have 
been  concerned  with  both  surface  and  body  waves.  We  have  achieved  further  success  in  predict¬ 
ing  the  refraction  and  multipathing  of  Rayleigh  waves  by  lateral  variations  in  crustal  structure. 
Using  body  waves,  we  have  completed  a  refined  version  of  the  amplitude-distance  curve,  and 
discuss  the  nature  of  precursors  to  pP,  S,  and  SKS  from  certain  deep  focus  events.  Two  other 


studies  describe  body  wave  scattering,  one  as  upper-mantle  structure  affects  PcP/P  amplitude 
ratio  's,  and  the  other  as  crustal  inhomogeneities  affect  the  P-wave  coda  at  short  distances. 

finally,  we  outline  three  investigations  that  are  continuations  of  earlier  work.  Our  study 
of  the  tectonics  of  Asia  includes  a  detailed  analysis  of  the  Baikal  Rift  Zone;  we  conclude  our 
analysis  of  NORSAR  travel  time  anomalies  by  discussing  their  extension  into  aseisinic  regions, 
and  we  show  an  application  of  the  single-channel  event  detector  described  in  the  preceding  3ATS. 

We  are  developing  a  sophisticated  interactive  computer  facility  as  an  in-house  research 
capability  that  will  enable  us  to  tccess  and  manipulate  seismic  data  c;  ce  they  become  available 
on  the  ARPANET.  Progress  in  the  design  of  the  appropriate  software  is  continuing,  and  will  be 
reported  in  more  detail  in  the  next  SATS. 

M.A.  Chinnery 
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SEISMIC  DISCRIMINATION 


1.  INTERNATIONAL  SEISMIC  MONTH 

A.  INTERNATIONAL  SEISMIC  MONTH  (ISM)  FINAL  EVENT  LIST 

Final  adjustments  have  been  made  to  the  ISM1'2  Event  List.  The  resulting  list  of  996  events 
and  some  discussion  of  its  statistics  have  been  published  in  a  recent  Technical  Note. 

The  ISM  detection  and  location  threshold  has  been  estimated  world-wide  and  for  Asia.  For 
this  purpose,  Asia  was  defined  using  standard  seismic  regions.4  Specifically,  we  took  Asia  to  be 
composed  of  the  seismic  regions  listed  in  Table  1-1.  Figure  1-1  shows  a  cumulative  histogram 
of  the  m.  values  for  the  859  ISM  events  for  which  we  were  able  to  determine  magnitudes  and  for 
the  subset  of  events  in  Asia.  In  the  Asian  region  there  were  215  events  with  mb  assigned,  and 
33  with  no  mb  determined. 


SEISMIC 

TABLE  1-1 

REGIONS  DEFINING  ASIA  FOR  THE  ISM 

Region 

Region  Name 

No. 

19 

Japan -Kuriles -Kamchatka 

26 

India -Tibet -Szechwan- Yunan 

27 

Southern  Sinkiang  to  Kansu 

28 

Alma-Ata  to  Lake  Baikal 

29 

Western  Asia 

30 

Middle  East-Crimea-Balkans 

41 

Eastern  Asia 

42 

N.  E.  Asia,  Northern  Alaska  to  Greenland 

47 

Baluchistan 

48 

Hindu  Kush  and  Pamir 

49 

Northern  Asia 

It  is  often  assumed  that  the  log  of  the  number  of  events  for  a  fixed  time  interval  in  any  region 
will  vary  linearly  with  mb>  One  can  then  fit  straight  lines  to  histograms  such  as  those  on  Fig.  1-1 


and  estimate  at  what  magnitude  level  the  detection  system  seems  to  start  missing  events.  Basham 
and  Anglin5  used  preliminary  ISM  data  for  a  region  similar  to  our  Asian  region  and  concluded 


that  the  cumulative  90-percent  detection  threshold  for  that  region  was  about  4.0.  Since  incremen¬ 
tal  thresholds  are  normally  0.2  to  0.3  units  higher,5  they  estimated  the  90-percent  incremental 


threshold  to  be  4,^  for  that  region.  The  data  of  Fig.  1-3  do  not  contradict  the  conclusion,  although 
this  fitting  of  straignt  lines  must  be  accepted  with  some  skepticism  as  a  method  of  estimating  De¬ 
tection  thresholds  since  it  does  use  a  theoretically  unsubstantiated  assumption  about  the  linear 
variation  of  seismic  activity  with  nr^. 

Fitting  the  complete  ISM  cumulative  occurrence  of  events  to  obtain  a  detection  and  location 
threshold  is  doubly  bothersome.  Not  only  must  the  straight-line  assumption  be  accepted,  but  it 
is  even  more  difficult  to  find  a  straight  region  to  fit  with  a  straight  line.  However,  if  such  a 


> ..  >-W. .  1. 


*m*w^ 


7***P^4, 


1  *^»^y  , ,  _ , 


method  is  used,  it  is  our  estimate  that  the  90 -percent  cumulative  threshold  obtained  is  between 
4.5  and  4.6,  with  the  corresponding  incremental  threshold  0.2  to  0.3  mb  units  higher. 

These  detection  capabilities  have  been  iurther  substantiated  by  comparisons  with  LASA  and 

NORSAR  bulletins.  Details  will  be  found  in  Lacoss,  et  al. 

Only  recently  has  the  Bulletin  of  the  International  Seismological  Center  (ISC)  been  published 
for  the  ISM  time  period.  Since  this  is  the  most  comprehensive  regular  bulletin  prepared  at  this 
time,  it  is  of  interest  to  compare  it  with  the  ISM  list. 


TABLE  1-2 

DISTRIBUTION  BY  GRADE  AND  mfc,  OF  THE 
WHICH  DID  NOT  APPEAR  ON  THE  ISC 

481  ISM  EVENTS 
BULLETIN 

Grade 

No.  of  Events 

™b 

No.  of  Events 

A 

49 

0.0 

84 

AI 

35 

<3.5 

44 

B 

17 

3.5-3.59 

25 

BI 

9 

3.6-3.69 

30 

C 

27 

3.7-3.79 

27 

Cl 

22 

3.8-3.89 

34 

D 

137 

3.9-3.99 

35 

DI 

185 

4.0-4.09 

39 

4.1-4.19 

32 

4.2-4.29 

32 

4.3-4.39 

28 

4.4-4.49 

26 

>4.5 

45 

1  3 

The  ISC  reported  515  of  the  996  ISM  events  (51.8  percent).  The  distribution  of  grade  ’  and 
mb  for  the  ISM  events  not  reported  by  the  ISC  is  given  in  Table  1-2.  Conversely,  the  ISC  Bul¬ 
letin  reported  299  events  during  this  period  that  did  not  appear  on  the  ISM  list.  Of  these  299 
events,  76  had  3  or  less  reporting  stations,  72  had  their  farthest  reporting  station  within  3  of 
the  epicenter,  and  113  events  did  not  have  more  than  3  of  the  stations  used  in  the  ISM  study  among 
their  reporting  stations..  A  breakdown  of  the  remaining  38  events,  of  which  3  events  (denoted 
by  *)  were  reported  on  the  ISM  unverified  array  event  list,3  that  possibly  could  have  met  ISM  ac¬ 
ceptance  criteria  is  shown  in  Table  1-3.  The  magnitudes  listed  in  this  table  are  of  various  kinds 
as  reported  by  the  ISC  with  only  the  teleseismic  body  wave  magnitude  mb  being  directly  com¬ 
parable  to  the  ISM  magnitude.  The  others  are:  ML,  local  magnitude;  M,  magnitude  (unspecifieci); 

and  M  surface-wave  magnitude.  In  general,  it  would  appear  that  no  events  above  the  claimed 
s 

ISM  detection  thresholds  appeared  on  the  ISC  Bulletin. 

There  are  two  primary  reasons  why  these  38  ISC  events  did  not  appear  in  the  ISM  event  list. 

In  some  cases,  arrivals  were  associated  but  we  rejected  the  resulting  hypocenter  because  of 
poor  residuals  or  insufficient  control  on  the  hypocenter.  More  commonly,  we  failed  to  ever  asso¬ 
ciate  the  arrivals  with  a  single  event. 
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TABLE  1-3 

ISC  BULLETIN  EVENTS  NOT  CONTAINED  IN  THE  ISM  LIST.  ALL 
ISC  EVENTS  WITH  3  OR  MORE  ISM  STATIONS  REPORTING  AND  AT 
ONE  STATION  MORE  THAN  3°  FROM  THE  EPICENTER  ARE  INCL 

MISSED 

LEAST 

UDED. 

Location 

No.  of  ISC 
Reporting 
Stations 

No.  of  ISM 
Stations 

Greatest  A 
of  ISC 
Stations 
(deg) 

Magnitude 

Austria 

13 

7 

3.96 

1.6  ML 

New  Hebrides 

4 

4 

4.80 

3"8  mb 

New  Hebrides 

6 

4 

5.62 

3.4  mb 

Off  Coast  No.  California 

8 

4 

9.32 

3.8  ML 

New  Hebrides 

5 

4 

4.10 

3.7  mb 

New  Hebrides 

5 

4 

6.84 

4.0  mb 

New  Hebrides 

4 

4 

4.63 

3.4  mb 

New  Hebrides 

5 

4 

7.60 

4.0  mb 

So.  Africa 

6 

4 

10.05 

3.3  M 

New  Hebrides 

9 

8 

151.51 

3.9  mb 

Norwegian  Sea 

7 

5 

10.97 

— 

So.  Africa 

6 

5 

9.41 

3.5  M 

So.  Africa 

8 

5 

11.00 

4.0  m, 
b 

E.  Caucasus 

13 

5 

65.56 

3.6  Ms 

New  Hebrides 

6 

5 

32.95 

3.4  mb 

Hindu  Kush 

6 

5 

74.72 

3.7^ 

No.  California 

7 

4 

5.54 

3.5  ML 

*  Greece 

11 

3 

25.61 

3.1  ML 

No  California 

9 

6 

21.27 

3.8  ML 

Hindu  Kush 

5 

4 

42.96 

4.2  rrih 

Off  Coast  No.  California 

9 

6 

12.59 

4.1  ML 

Davis  Strait 

5 

4 

35.03 

Western  Caucasus 

10 

4 

26.57 

3.7  M 

s 

No.  California 

11 

6 

10.39 

3.5  ML 

*  Hindu  Kush 

7 

6 

44.52 

3.7  mb 

New  Hebrides 

8 

7 

148.05 

4.0  m, 

b  | 

Off  E.  Coast  New  Zealand 

15 

6 

151.17 

4.4  M 

No.  Island  New  Zealand 

22 

7 

146.36 

5.0  M 

Loyalty  Isles  Reg. 

11 

8 

98.66 

4.7 

*  Chile -Argentina  Border  Reg. 

7 

6 

83.88 

— 

Off  Coast  No.  California 

8 

4 

7.81 

4.1  ML 

California-Mexico  Border  Reg. 

9 

5 

19.34 

4.5  ML 

Off  Coast  No.  California 

7 

4 

9.90 

3.9  ML 

So.  Africa 

15 

11 

146.11 

4.1  M 

Loyalty  Isles  Reg. 

5 

4 

67.38 

3.8  n^ 

Loyalty  Isles  Reg. 

5 

4 

33.37 

3.8  mb 

New  Britain  Reg. 

10 

9 

25.06 

4.9  ML 

Tadzhik 

5 

4 

42.96 

3 
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As  explained  earlier,  the  location  programs  used  for  the  ISM  find  the  origin  time  and 
hypocenter  which  minimize  a  weighted  sum  of  the  squares  of  observation  residuals.  An  error 
ellipsoid  is  defined  which  corresponds  to  a  unit  increase  of  this  sum  of  squares.  From  this 
ellipsoid,  we  obtain  two  measures  of  epicenter  quality.  One,  MAXAX2  is  the  major  semiaxis 
of  the  epicenter  error  ellipsoid.  The  other,  DEPTHQ,  is  half  the  vertical  extent  of  the  hypo- 
center  error  ellipsoid. 

Experiments  have  been  completed  to  calibrate  these  quality  variables.  In  one  of  these  ex¬ 
periments,  all  well-recorded  events  which  used  depth  phases  in  obtaining  a  hypocenter  were  re¬ 
located  without  using  the  depth  phases.  Of  these,  141  converged  properly  without  any  need  for 
restraining  depth  for  the  relocation.  Since  the  previous  .  (cations  with  depth  phases  tended  to 
pin  the  depth  closely  to  that  indicated  by  the  depth  phases,  we  took  the  difference  between  the 
two  depths  as  the  actual  error  in  depth  for  the  relocated  events.  The  absolute  value  of  this  error 
was  weighted  inversely  by  the  depth  quality  of  the  relocated  event,  and  the  histogram  of  Fig.  1-2 
was  generated.  Accepting  the  assumption  that  the  original  depth  was  correct,  we  see  that  the 
indicated  depth  of  the  relocated  events  is  almost  never  in  error  by  more  than  2*  DEPTHQ  kilom¬ 
eters  and  that  it  is  less  than  DEPTHQ  for  75  percent  of  the  relocated  events.  Other  calibration 
experiments,  involving  MAXAX2  as  well  as  DEPTHQ,  are  described  in  Ref.  3.  The  results  are 
all  similar  to  those  indicated  on  Fig.  1-2.  R  T  Lacoss 


R.  E.  Needham 


B.  PERFORMANCE  OF  SELECTED  SHORT-PERIOD  SITES 
DURING  THE  ISM  PERIOD 


With  the  completion  of  the  ISM  experiment,  there  is  now  available  a  large  quantity  of  station 
data,  both  short-  and  long-period  arrival  times,  amplitudes,  etc.  A  great  deal  of  effort  has  gone 
into  merging  the  short-period  data  for  over  30  selected  sites  with  the  epicenter  location  data  from 


the  ISM  list.  The  total  quantity  of  data  exceeds  350,000  items,  all  contained  on  one  magnetic 

8 

tape  in  a  format  acceptable  to  the  DADS.  With  the  data  in  this  format,  it  is  now  possible  to  dis¬ 


play  all  or  part  of  it  with  the  DADS.  In  addition,  the  complete  library  of  routines,  both  numerical 
and  graphic,  are  available  to  aid  in  the  analysis  of  single  station  or  sets  of  stations  or  in  epicenter 
location  or  any  other  related  phenomenon. 

The  analysis  of  some  of  the  station  performance  characteristics,  which  is  just  beginning, 
will  cover  all  the  stations  that  reported  a  substantial  number  of  arrivals  during  the  ISM  period, 
with  an  emphasis  on  the  Canadian  stations.  Initially,  only  the  better  events,  i.e.,  those  reported 
with  a  quality  of  A  or  AI,  will  be  used  in  the  study.  This  will  give  a  more  reliable  base  of  events 
that  have  a  large  number  of  stations  used  for  locations  and  magnitude  determinations.  For  these 
375  events  with  A  or  AI  rating,  the  station  detection  levels  are  shown  in  Fig.  1-3  where  it  can  be 
seen  that  only  7  stations  detected  more  than  half  of  the  375  events.  Of  these  7  stations,  4  are 
actually  seismic  arrays  whose  detection  levels  were  greatly  enhanced  by  beamforming.  The  re¬ 
maining  3  stations  (MBC,  UBO,  and  KBL)  are  simply  good  quality  single  site  seismometers. 

The  Canadian  seismic  network  data  that  contributed  to  the  ISM  provided  a  very  good  quality 
set  of  data  from  a  single  area.  This  will  provide  a  good  data  set  to  examine  variables  such  as 
travel  time  anomalies,  magnitude  correction  terms,  and  other  related  factors.  One  term,  the 
magnitude  correction  term  in  the  formula  for  short-period  magnitude,  has  been  studied  initially. 
While  there  are  many  ways  to  determine  what  the  magnitude  correction  values  are,  the  method 
used  here  was  to  plot  magnitude  of  all  A  and  AI  events  for  a  given  station  on  a  distance/magnitude 


. . —  . 


plot.  Then  the  outline  of  the  Gutenberg-Richter  magnitude  "B"  factor  curve  for  depth  0  and  700 
was  overlayed  onto  the  data  and  fitted  to  the  data  by  eye.  The  resulting  value  on  the  0  depth 
curve  at  55°  was  chosen  as  the  magnitude  detection  level  for  that  station.  The  results  for  the 
Canadian  network  are  displayed  on  Fig.  1-4  along  with  possible  contours  of  the  magnitude  detection 
levels  throughout  all  of  Canada  based  on  the  network  values.  Only  one  site,  SES,  shows  a  radical 
departure  from  the  pattern.  It  has  an  unusually  low  magnitude  detection  level.  There  can  be 
several  reasons  for  this  but  it  is  probably  due  to  the  very  low  detection  level  to  the  Gulf  of 
California  swarm  in  the  ISM. 

Continued  examination  of  all  factors  relating  to  the  ISM  list  and  the  station  set  hopefully  will 
reveal  more  interesting  features  in  the  data  and  possible  explanations  of  these  effects. 

R.  M.  Sheppard 
R.  E.  Needham 

C.  LONG-PERIOD  RESULTS  FROM  THE  INTERNATIONAL  SEISMIC  MONTH 

The  first  study  of  the  long-period  aspects  of  the  ISM  has  been  completed  and  reported  else- 
g 

where  ;  here,  we  shall  only  review  some  of  the  problems  encountered  and  the  results  of  the  study. 
The  chief  problems  were  found  in  the  automatic  association  of  the  long-period  detections  with  the 
epicenter  list  and  in  the  surface  wave  magnitude  assignments.  Only  half  of  the  detections  were 
assigned  an  association  by  the  individual  readers,  persons  with  widely  varying  experience  in  the 
interpretation  of  long-period  seismograms.  Therefore,  to  treat  all  the  data  in  a  consistent  man¬ 
ner,  it  was  decided  to  consider  all  the  long-period  detections  as  unassociated  and  develop  various 
automatic  association  schemes.  The  definition  of  what  is  an  optimum  scheme  is  somewhat  arbi¬ 
trary;  a  scheme  that  makes  no  misassociation  will  probably  give  too  few  true  associations  to  be 
practical,  while  another  that  associates  all  the  readings  may  make  an  unacceptable  number  of 
mistakes.  Various  combinations  of  path  dependent  group  velocity  curves  and  association  win¬ 
dow  lengths  were  tested;  however,  it  was  found  that  a  world  average  group  velocity  curve  with 
a  ±0.5 -km/sec  association  window  gave  more  unique  associations  than  any  other  scheme.  This 
latter  scheme  also  gave  the  best  results  when  compared  with  the  cases  where  an  association  was 
made  by  the  readers.  In  all,  about  55  percent  of  all  the  detections  identified  as  main  Rayleigh 
phases  were  given  unique  associations  by  the  accepted  scheme. 

The  "Prague"  formula  for  surface  wave  magnitude  determination  was  used  initiaUy,  but  it 

was  found  to  give  consistent  results  only  near  periods  of  20  sec,  where  it  was  defined.  The 

10 

method  described  by  Marshall  and  Basham,  which  uses  a  path  dependent  period  connection,  was 

then  applied  and  gave  much  better  results,  i.e.,  it  yielded  consistent  magnitudes  for  an  event 

over  a  wide  range  of  periods  and  distances,  in  order  to  incorporate  the  maximum  number  of 

long-period  detections  into  the  magnitude  determinations,  the  distance  tables  of  Marshall  and 
10 

Basham  were  extended  to  180°  and  the  period  corrections  were  extended  from  40  to  100  sec. 

By  using  these  extended  tables,  individual  Mg  values  were  computed  for  each  associated  long- 

period  detection,  and  an  event  average  Mg  and  standard  deviation  were  computed. 

Figure  1-5  shows  the  distribution  of  these  average  M  values  plotted  as  a  function  of  the  m, 

S  D 

of  the  grade  A  events.  This  grade  was  used  to  indicate  the  more  accurately  located  events  of  the 
ISM  short-period  list.  The  curve  represents  the  distribution  of  world-wide  M  -m,  data  suinma- 

11  I?  8  “ 

rized  by  Aki,  and  the  straight  line,  due  to  Gutenberg  and  Richter,  is  often  used  as  an  M  -m, 

S  D 

relationship.  It  is  evident  that  this  line  is  not  representative  of  either  Aki's  summary  curve  or 


the  world-wide  ISM  data.  Figure  1-6  gives  some  indication  of  capabilib  of  the  l 

W°rk  “Sed  “  ,he  BM  ‘°  «=*e=t  surface  waves  from  shallow  events  In  Asia 

number  of  events  of  a  tri^n  ™  r  events  m  Asia.  It  shows  the  total 

magnitude  based  on  two  or  mo^m^a^  ^  ^  ^  nUmber  °f  that  total  for  which  a  surface  wave 

percent  cumulative  threshold  for  M  ^  °"  ^  ^  ^  9°' 

events  during  the  ISM  was  about  m  S=  4  6  F  measurements*  for  shallow  Asian 

b  ‘  ■ 

J.  R.  Filson 
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Fig.  1-1.  Cumulative  histogram  of  complete  ISM  mb  and  for  events  limited 
to  Asian  seismic  regions  19,  26,  27,  26,  29,  30,  41,  42,  47,  48,  49. 
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normalized  change  in  depth 

Fig.  1-2.  Histogram  of  change  in  depth  resulting  from  relocating  A  grade 
events  after  discarding  all  depth  phases.  Each  change  is  normalized  by  re¬ 
location  depth  quality  parameter.  Total  number  of  events  used  is  141. 
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Fig.  1-3.  ISM  event  detections  at  selected  stations. 


II.  DATA  MANAGEMENT 


GLOBAL  SEISMIC’  DATA  MANAGEMENT  SYSTEM 

Dro.eT  tAhRPA  NMROhaS  data  collection  and  utilisation 

project  This  project  tnvolves  the  deployment  of  addition,,  seismic  instrumentation-  acquisi- 

ton  of  digit  at  data  real  time  and  via  tape  recordings;  ,he  ,„alysis  of  such  data  to  obtain  event 

Hat,  and  o  tes,  dtscrin, inants  and  operational  methods;  and  the  organisation,  ,torage  retrieva, 
an  use  of  raw  and  processed  data  in  research  programs  directed  toward  improving  the  „  S  ' 
seistmc  discrimination  capabilities.  In  the  sequel,  we  will  refer  to  the  system  which  is  gomg 
ereaise  as  the  Integrated  World-Wide  Seismic  System  (1WWSS).  The  lesign  and  imple- 
men  , uon  of  1WWSS  is  a  complex  program  requiring  careful  coordination  to  achieve  a  tin.  ly 
and  effective  working  syst.  m  whose  characteristics  have  been  optimised  toward  the  solution  of 
remaining  seismic  problems  related  to  discrimination.  To  this  end.  Lincoln  has  started 
working  close, y  with  NMlto,  other  Oovernmen,  agencies,  and  contractors  be  certain  tha,  the 

.To"'!",  Se,Sn  <CSlly  “SefUl  *to*  *«  <•  eood  coordination  and  communica! 

tion  between  participants  and  subsystems. 

We  discussed  th  global  seismic  data  project  at  length  with  NMRO  and  participants  in  the 
projec  to  obtain  from  each  his  concept  o,  the  system  and  re.ated  information  Emphlsis  was 
n  e  kind  and  quantity  of  seismic  data  and  the  overall  flow  and  management  of  the  data  An 

LTaTnTiut  ;:rem  which  is  evoivine  was  prepared  and  discussed  ^  ^0^  t0  be 

tain  n.it  there  were  no  major  discrepancies  or  misunderstandings 
One  niaior  component  of  the  IWWSS  will  be  a  mass  store  device  for  processed  and  unproc- 

,o  he  r'T  dT’  A  WOrk‘"g  d°CUmCnl  ^  PrCPared  ”hlch  characterizes  the  major  files 

da,  n  ‘  r<!la,i0nships  *h«e  files  and  other  1WWSS  component,  and  the 

ow  into  and  between  these  files.  This  document  and  considerable  other  information  has 
n  entered  into  the  On-Line  System  which  is  discussed  below. 

Lincoln  has  started  ,o  use  the  On-Line  System,  NLS,  developed  by  Stanford  Research  In 
;  1,  u  e.  as  a  working  fool  for  the  eons, rue, ion  of  a  directory  that  will  eventually  1,^1  d- 
rable  informal, on  about  the  1WWSS.  The  NLS  is  supported  by  the  TEN11X  system  at  Office-1 
on  the  AREA  network.  The  directory  llself  win  contain  many  files,  each  of  Lch  wm  con¬ 
cerned  w„h  specific  , epics  of  ,he  integrated  world-wide  system.  One  file  OVERVIEW  „  , 
contain  a  general  description  of  the  system  and  system  components.  Within  that  file  there  will 
ink  references  (an  NLS  feature)  to  more  specific  documents  which  will  contain  more  detailed 
informauon.  There  will  probably  be  links  throughout  all  documents  to  specific  and  more  detailed 
descriptions  of  the  areas  of  interest.  When  the  directory  is  completed,  an  ARPANET  user  will 
be  able  to  view  a  general  description  of  the  IWWSS  and  then  be  able  to  jump  to  more  detailed  de¬ 
scriptions  and  continue  to  quickly  cross  -examine,  thn  •  f 

of  experience  with  .K.tn  informal, on  available.  Only  a  small  amount 

umel  Ze7  7  Perf°rm  *  am, nation  of  the  availabfe  doc- 

users  wit  a  Zd  T7  ““  dlreC‘°ry  W‘“  Pr°Vide  SyS‘em  d»'lpp^-  and  evenlually 
users,  with  a  raptd  access  to  descriptions  of  current  concepts,  processors,  software,  hardware 

01  ganizations,  and  people  involved  in  the  IWWSS.  ’  ' 

severalT1™  SkeIet°"  "  OVERVIEW  "'arly  completed,  along  with  the  definition  of 

sever  ,  other  supporting  files.  Generally,  ,he  first  NLS  stalemen,  in  each  file  will  be  a  .able 

Of  contents  o,  tha,  file,  wh.le  ,he  second  NLS  statement  will  contain  all  the  links  that  wer  c Red 
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Within  the  file.  Vhe  table  of  contents  is  provided  so  that  general  references  to  the  file  via  link, 
from  other  docum-nts  can  supply  the  user  with  a  quick  look  at  the  contents  of  the  file.  The  link 
re  e  -cnees  cited  within  a  file  win  be  contained  in  ,he  second  statement  of  that  file  so  tha,  if  ere 

^“1  ZTrTl  al,era“0nS  *  C‘ted  d°C“™'-  *h“  “  -»  -  ~r  to  find  and  eh  J  I 

mLs  store  wh  h  d”  “  al'"'d  d°CU"Cnl'  cteoac.eriai„g  major  files  on  the 

,  ic  wis  iscussed  above,  is  currently  part  of  this  OVERVIEW  Tie 

varie.yrof  LwlT  TT  ‘  ^  •>*“  tes  >“0”  designed  to  facilitate  a 

outside  o  ^  “P  “  "»  I"—  «">««■  application 

d,°f  h:  """***“  oommunity  of  developers.  Nevertheless,  we  .Clcipate  tha,  i,  will  be 

1WWSS  “  “  J  '°r  lhe  dc''el°i,m'!nt.  management,  and,  eventually,  final  documentation  of 
IdThe  ,  ’  *he  SelSm‘C  "Ser  sl»-..y  will  need  to  be  expanded  Pe¬ 
so  that  we  ma  ”  We  are  Cl!rrently  organizing  the  seismic  information  within  NTS 

so  that  we  may  expand  the  user  group  and  establish  closer  dialogs  and  coordinate  the  end 

1  “I! :tr“C,Ur'  “f  ,hC  ““  id“,ly  —  *  tnntnal  activities  by  van-’ 

R.  T.  Lacoss 
R.  M.  Sheppard 
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in.  earthquake  source  mechanisms 


A.  FREQUENCY -MAGNITUDE  CURVES  AND  THE  RELATIONSHIP 

appuiir ~~  rr  vai“bi° -  *■ 

::  «„e  p.otS  r„  fairly  ,he  ^ 

NOAA  ’  a„mdbUSGS,hrV0,Heen  ° ^  ^  ° f  .produced  bv  OSCGS 

.ho^b  ,ei"'S  10  i9”-  Subs=*a  «'  Mo  lata  set  show  similar  features  o 

or  „,b  ;r:,  •*  tis  °r  -  —  -  ***,«. 

assumed.  If  a  straight  line  is  passed  through','  he'data  ZZc  <  rT<T  '"h  “  “ 

a  b-va.ue  0,  1.35,  which  is  very  high.  b„,  agrees  with  the  results  of  Everndem'  *  ‘  “  “ 

suiting  'curvTis  PDF  «u  T  ^  ^  °'  '"V  The  bo¬ 
at  n,  <  5  ,f  ,1  ntagnitudes  greater  Utan  5.  bu,  shows  a  divergence 

ity)  bthc  points  ma  ‘  7brVCS  SrC  matChed  f°r  mb  *  5  (t°  allOW  f°r  Slight  variations  in  tcial  activ- 

«™l,h.  Z  .*  ,  y  *  "OSS  °"  lri6' IIM  Sh»  »*>  -ata.  The  trend  indicated  by  the 

Plate  at  n,  -  4  5  Z  O  Ca,‘u'u‘  hck,w  mb  '  5’  s‘”«  ««'  ISM  data  were  certainly  not  com- 

b  •  ,  i&  trend  probably  continues  at  least  down  to  m  -  4  5  Th„ 

tween  the  two  data  sets  is  presumably  the  result  incompleteness^  the  PDF  iTT™* 
below  m,  r  5.  Ther»  is  no  uviHAncn  th  .  „  ieteness  in  the  PDE  data  at  magnitudes 

as  has  been  ““  ^  ""  *°  *  ■"»“»  »—■  *■  *-  «*■. 

Because  M  is  not  assigned  to  all  events  in  the  pnp  u-«  •.  • 

I  be  b-value  for  the  frcquency-M  curve  is  0.92,  in  good  agreement  with  *  4. 

v  id  e  ^c  on  sis  t  enT  rn  e  f  ^  °"  ^  ^  “  aPPearS  that  "l 

a  consistent  measure  of  seismic  moment  for  magnitudes  up  to  about  7S  At  L‘  ,  s 

the  effect  of  corner  frequency  is  felt,  and  Gutenberg  and  Richter-s  data  suggest  thaTthe IroZe  ’  - 

M  curve  becomes  vertical  in  the  vicinity  of  M  -  8  5  M  v  1  n  a  qu  ncy 

should  setdom  exceed  this  value.  regardLl  ^seisnde  “iu^ ZT.  * 

ment  with  source  mechanism  theory.3  g  nera  agrec- 

fhe  slope  of  the  m,  curve  indicatos  that  m  1=  „  * 
ror  the  range  o,  the  data*  shown,  nr  »  4  “ 1 ^  T”"  *  SeiS'"‘C  m°mM 

tarn  over  a  slope  0,  about  0.,  it  ls  no,  ‘T 

~  r  “■  <™»  ™b  >  4.  A,  its  upper  end.  the  frequel  "  '  a~ 

off.  As  in  the  ease  of  the  Ms  data,  this  is  due  to  a  corner  frequency  effect  and  ,H 

to  become  vertical  at  m  ~  6  7  m  ,  1  i  q  y  ffect>  and  the  curve  appears 

questionabte.  b  b  *  “  larg°r  'ha"  ^  ^  llk'h  to  be  ex, rente, y 

Althougl  there  arc  some  problems  in  attempting  compare  the  two  curves  in  Fig  111  ,  be 

cause  they  were  drawn  from  different  sets  of  earthquakes  it  is  nncCM  *  g‘  1  ’  be_ 

the  general  form  of  the  world  average  M  -m  relation  Cle  ^  1  Gm  t0  diSCUSS 

4  5  <  m  <  k  .*u  ,  ,  •  S  s  b  relatl0rl-  Clearly,  this  will  be  a  straight  line  for 

b  -  ,  wi  a  slope  that  is  easily  calculated  to  be  1.43.  At  its  upper  end  for  m  >6  the 

corner  frequency  effect  causes  the  Ms-xnb  curve  to  rise  toward  the  vertical.  ’  5  ’ 
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An  Mg-m^  relation  derived  directly  from  the  two  curves  in  Fig.  Ill— 1  is  shown  in  Fig.lII-2. 

I  ossible  inferences  in  the  definition  of  Mg  could  move  this  curve  up  or  down.  Note,  however, 
that  it  shows  good  agreement  with  the  overall  scatter  in  Mg-mb  observations  quoted  by  Aki,3  and 
with  the  ISM  data  given  by  Filson.4  Slightly  better  agreement  occurs  if  the  curve  is  raised  by 
0.2  unit  of  Mg,  and  this  is  eonfirned  by  what  PDF  data  are  available  (data  not  shown). 

’’’he  Gutenberg-Richter  relation 

Ms  =  1.59  mb  -  3.97 

is  a  very  poor  fit  to  the  data,  and  should  be  abandoned.  A  reasonable  fit  to  the  curve  in  Fig.  Ill -2, 
moved  vertically  up  by  0.2  unit,  is  given  by 

Ms  =  2.91  -  0.68  mb  f  0.215  mb  for  4. 5  <  <  6.5  . 

1  ntil  more  Mg  data  are  available,  this  relation  should  be  considered  as  a  useful  approximation 
to  the  world  average  Mg-mb  relationship.  Significant  departures  from  this  expression  are  to 
be  expected  within  limited  geographical  regions. 

M.A.  Chinnery 
R.  G.  North 

11.  TIIE  MOMENT-Ms  RELATIONSHIP,  AND  THE  FREQUENCY 
OF  LARGE  EARTHQUAKES 

Estimates  of  seismic  moment  are  now  available  for  a  large  number  of  earthquakes.  We 
have  selected  these  data  for  events  with  a  reliable  value  of  Mg.  These  are  either  published 
PDF.  Mg  values,  where  a  considerable  number  of  stations  were  used  in  the  determination,  or 
Mg  values  determined  directly  from  the  spectral  amplitudes  of  surface  waves  as  the  moment 
was  determined.  A  compilation  of  91  MQ-Mg  pairs  is  shown  in  Fig.  111-3,  and  includes  ail  re¬ 
liable  events  for  M  >5. 

s 

The  shape  of  the  MQ-Mg  relation  provides  another  valuable  constraint  for  source  mechanism 
theories,  and  it  is  therefore  of  interest  to  attempt  to  use  the  data  points  to  define  a  mean  rela¬ 
tionship.  First  we  note  that  the  data  define  a  remarkably^  straight  line,  when  the  three  largest 
moment  values  are  omitted.  This  line  has  the  equation 

log  Mq  =  16.29  +  1.5  Mg  (5<  Mg  <  8) 

This  line  is  a  considerably  better  fit  to  the  data  than  the  w2  model  of  Aki3  which  is  also  shown 
in  Fig.  Ill  -  3 . 

The  upper  end  of  this  curve  is  of  considerable  interest,  although  it  is  only  poorly  defined 
by  the  present  data.  W'e  expect  a  corner  frequency  effect  to  dominate  at  some  magnitude,  and 
the  Mq-Ms  curve  to  trend  toward  the  vertical.  The  data  suggest  the  presence  of  this  trend,  but 
do  not  define  it  elearly. 

One  unusual  approach  to  this  data  set  is  suggested  by  the  frequency-M  eurve  shown  in 

s 

Fig.  111-4.  Can  the  departure  of  the  frequency-Mg  curve  from  a  straight  line  (at  M  >  8)  be  at¬ 
tributed  entirely  to  a  corner  frequency  effect?  If  so,  this  would  imply  tha„  the  frequency-moment 
relation  should  be  a  straight  line.  Such  a  suggestion  has  been  made  by  Wyss,5  although  he  gives 
no  justification  for  the  hypothesis. 

We  can  reverse  the  argument  as  follows:  Suppose  that  the  frequency-moment  graph  is  a 
straight  line,  what  is  the  implied  moment-Mg  relationship?  Clearly,  the  moment-Mg  relation 
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must  be  a  straight  line  up  to  Mg~  8,  and  it  will  then  curve  upward.  The  solid  line  in  Fig.  I1I-3 
shows  such  a  relationship,  and  the  corresponding  frequency-moment  relation  is  shown  in 
Fig.  Ill— 4. 

Although  some  uncertainties  in  the  Mg  values  used  in  Gutenberg  and  Richter's  data  remain, 
we  appear  to  be  able  to  draw  the  following  conclusions: 

(1)  The  available  data  are  entirely  consistent  with  the  hypothesis  that  the  frequency- 

moment  relation  is  a  straight  line  out  to,  and  probably  beyond,  M  =  10^  dyne-cm 

and  that  the  efiect  of  corner  frequency  starts  to  affect  Me  determinations  (20-sec 

s 


period)  at  just  about  Mg  =  8.0. 
abie  for  any  earthquake. 


M 


values  of  greater  than  8.5  are  probably  unreli- 


(2) 


,  for  M  >  8,  has  been 

O  O 


c. 


The  fall -off  in  frequency  of  large  earthquakes  defined  by  M 

used  as  an  argument  that  earthquakes  with  source  dimensions  larger  than  those 
con  monly  attributed  to  Mg  =  8.5  are  impossible  (see,  for  example,  the  arguments 
of  Ilousner  ).  The  frequency-moment  curve  in  Fig.  111-4  shows  clearly  that  there 
is  no  foundation  for  this  argument.  It  is  possible  that  extremely  large  earthquakes 
<Mo  >  10  dyne-cm)  occur  from  time  to  time,  particularly  on  a  geological  time 
scale.  Although  there  may  be  a  maximum  possible  moment,  it  is  apparently  be¬ 
yond  the  range  of  presently  available  data.  The  whole  concept  of  the  "maximum 
possible  earthquake,"  as  commonly  used  in  earthquake  risk  studies,  is  very 
questionable. 

M.A.  Chinnery 
R.  G.  North 

DISTORTION  OF  APPARENT  EARTHQUAKE  FOCAL  MECHANISM 

BY  MANTLE  STRUCTURE  BENEATH  OCEAN  RIDGES 

Normal-faulting  earthquakes  on  the  crests  of  midocean  ridges  show  an  unusual  characteristic 
which  has  been  noted  only  in  passing  in  the  literature:  Well-constrained  fault-plane  solutions  de¬ 
rived  from  P-wave  first  motions  indicate  an  apparent  nonorthogonality  of  the  nodal  planes  when 
standard  projections  are  used  to  map  observations  back  to  the  focal  hemisphere.  A  particularly 
clear  example  of  this  phenomenon  for  first-motion  data  for  a  ridge-crest  earthquake  is  shown 
in  Fig.  lll-5(a-b).  This  earthquake,  which  occurred  on  20  September  1969  on  the  Reykjanes  ridge, 
is  obviously  a  normal-faulting  event.  The  angle  between  the  nodal  planes  in  the  diiatational  "quad¬ 
rant"  is  only  about  60°,  however.  Two  orthogonal  fault  planes  could  not  be  fit  to  the  data  without 
violating  a  large  number  of  them. 

It  is  our  contention  that  the  nonorthogonality  of  the  nodal  planes  is  an  artifact  of  the  projec¬ 
tion  used  to  map  the  earth's  surface  back  to  the  focal  sphere.  Specifically,  the  standard  tables 
of  ray  parameter,  or  angle  of  incidence,  vs  epicentral  distance  are  inadequate  when  the  source 
region  is  laterally  neterogeneous.  We  cannot  conclusively  rule  out  the  possibility  that  the  non- 
orthogonality  is  real,  e.g.,  that  the  earthquake  source  is  best  represented  by  a  double  couple  and 
a  superposed  explosive  component.  An  explosive  component  with  P-wave  amplitude  equal  to 
0.3  to  0.5  times  the  maximum  P-wave  amplitude  from  the  double  couple  would  yield  the  observed 
angle  of  60°  to  70°  between  the  best  fitting  nodal  planes.  In  particular,  a  source  consisting  of 
a  double  couple  and  a  superposed  explosive  component  with  amplitude  about  half  the  maximum 
doubl>  -couple  amplitude  is  entirely  consistent  with  the  first-motion  data  in  Fig.  1II-6.  Because 
of  the  clear  evidence  for  faulting  in  the  rugged,  blocky  topography  of  the  mid-Atlantic  ridge, 


however,  we  prefer  a  simple  dislocation  model  for  the  earthquake  mechanism  and  we  attribute 
the  observed  nonorthogonality  of  nodal  planes  to  a  path  effect. 

We  have  thus  constructed  models  of  the  upper  mantle  beneath  spreading  ridges  based  upon 
theoretically  calculated  temperature  fields  and  experimentally  determined  phase  relationships 
for  peridotite,  and  have  used  three-dimensional  ray-tracing  techniques  to  calculate  the  paths  of 
P-waves  through  these  models.  Figure  111-6  shows  examples  of  such  ray  paths,  which  can  be 
seen  to  be  bent  downward  by  a  region  of  very  low  seismic  velocity  beneath  the  ridge  axis  asso¬ 
ciated  with  extensive  melting  at  temperatures  above  the  dry  solidus  of  peridotite.  The  results 
of  these  calculations  can  easily  explain  the  observed  anomalous  focal  mechanisms,  as  shown  by 
the  corrected  data  of  Fig.  Ill- 5(b) . 

This  work  is  currently  being  extended  to  include  data  on  shear  wave  polarization,  which  we 

hope  will  be  able  to  conclusively  differentiate  between  source  effects  and  path  effects. 

B.  R.  Julian 
S.  Solomon 

D.  SOURCE  DIMENSIONS  OF  EVENTS  IN  WESTERN  TURKEY 

A  detailed  study  has  been  made  of  the  seismicity  of  Western  Turkey  and  the  mechanisms  of 
some  of  the  larger  events  in  this  area,  here  defined  as  that  part  of  the  country  west  of  32  °E, 
determined.  It  has  been  suggested  (McKenzie7)  that  the  extension?!  nature  of  the  deformation 
here,  geologically  well  evidenced  by  a  series  of  horst  and  graben  structures,  results  from  the 
motion  of  an  Aegean  plate  southwestward  away  from  a  Turkish  plate,  itself  moving  rapidly  west¬ 
ward  with  respect  to  Eurasia.  In  recent  years,  this  region  has  experienced  a  number  of  large 
and  destructive  earthquakes,  the  greatest  of  which  was  that  of  Gediz  on  28  March  1970. 

All  events  reported  in  this  area  for  1964-70  (ISC  bulletin)  and  1971-72  (Earthquake  Data  Re¬ 
ports  of  the  USGS)  have  been  relocated  using  a  program  by  Julian.  The  larger,  well-located 
events  were  used  to  define  a  set  of  station  corrections  which,  apart  from  the  closer  stations, 
appear  to  vary  little  across  the  area  considered.  These  corrections  were  then  used  in  a  "joint 
epicentre  determination"  method9  to  obtain  relative  locations  of  the  smaller  events.  The  re¬ 
located  epicenters  are  shown  in  Fig.  111-7.  The  aftershock  sequences  of  the  1967  (Mudurnu), 

1969,  1970  (Gediz),  and  1971  earthquakes  are  clearly  seen. 

A  study  of  the  aftershocks  for  the  1967  event  reveals  clearly  a  westward  migration  of  the 
shocks  with  time.  This  is  in  agreement  with  the  westward  progression  of  larger  events  along 
the  North  Anatolian  Fault  (at  the  westernmost  end  of  which  this  event  occurred)  first  noticed  by 
Ketin.10  The  aftershock  sequence  for  the  Gediz  event,  the  largest  in  this  area  in  the  last  five 
years,  is  shown  in  Fig.  II1-8.  The  trend  of  the  aftershock  distribution  agrees  well  with  that  of 
both  the  fault-plane  solution  determined  by  McKenzie7  and  the  surface  faulting  (Ambraseys  and 
Tchalenko11).  Its  length  is  somewhat  greater  than  that  of  the  latter,  and  the  distribution  shown 
yields  an  area  of  faulting  of  surface  dimensions  80  X  40  km.  Since  both  the  nodal  planes  of  the 
fault-plane  solution  dip  at  -45°,  the  actual  area  is  closer  to  80  X  60  km. 

For  the  events  of  1967,  1969,  and  1971,  the  aftershocks  define  the  area  of  faulting  much 
less  clearly.  Tentative  estimates  are  70  X  30,  50  X  20,  and  40  X  20  km,  respectively.  The  first 
agrees  well  with  a  length  of  surface  faulting  of  80  km  (Ambraseys  and  Zatopek  ). 

The  seismic  moments  of  10  events  in  this  region  prior  to  1971  have  been  determined  from 
the  amplitude  spectra  of  Rayleigh  waves.13  The  fault-plane  solutions  for  the  events  of  1971-72 
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are  virtually  identical,  viz.  normal  faulting  striking  100°  to  120°E  of  N,  with  those  of  1909-70 
(McKenzie,  personal  communication).  The  seismic  moments  for  the  former  can  thus  be  quickly 
estimated  by  measuring  maximum  amplitudes  of  Rayleigh  waves  at  a  number  of  statiom  "nd  com¬ 
paring  these  amplitudes  with  those  for  an  event  for  which  the  moment  is  already  well  deiei  uned, 
in  this  case  that  of  28  March  1970, 

The  source  parameters  for  16  events  in  Western  Turkey  are  given  in  Table  III— 1 .  The 

values  quoted  are  those  of  the  USCGS,  and  the  M  values  are  from  various  sources.  The  seismic 

s  24 

moment  IV1  and  its  standard  deviation  are  given  in  units  of  10  dyne-cm,  together  with  the  num- 
o 

ber  of  station  records  N  used  in  its  determination.  For  these  events  for  which  it  has  been  pos¬ 
sible  to  estimate  the  length  L  and  width  W  of  faulting,  we  may  estimate  the  mean  relative  dis¬ 
placement  u  taking  place,  through  the  definition  of  seismic  moment  M  =  puA,  where  p  is  taken 

11  ,2  ° 

as  3  X  10  dyne/cm  ,  and  A  is  the  area  of  the  fault  plane. 


TABLE  III-l 

EVENTS  IN  WESTERN  TURKEY  FOR  WHICH  THE 
MOMENT  Mq  HAS  BEEN  DETERMINED 

SEISMIC 

Date 

Latitude 

(°N) 

Longitude 

(°E) 

mb 

M 

s 

M 

0 

N 

L 

(km) 

W 

(km) 

u 

(cm) 

6/10/64 

40.40 

28.87 

5.9 

6.7  MOS 

180  ±  70 

13 

13/6/65 

37.89 

29.74 

5.1 

5.5  MOS 

8.2  ±  2.6 

6 

22/7/67 

40.73 

30.86 

6.0 

7.1  CGS 

1500  ±  860 

14 

70 

30 

80 

5/12/68 

36.63 

37.07 

5.4 

6.1  ATH 

18  ±  7.5 

9 

23/3/69 

39.24 

28.67 

5.6 

6.0  MOS 

9.1  ±  4.9 

15 

25/3/69 

39.16 

28.52 

5.5 

6.0  MOS 

19  ±  9.2 

19 

28/3/69 

38.86 

28.56 

5.9 

6.5  MOS 

120  ±  67 

16 

50 

20 

40 

28/3/70 

39.20 

29.57 

6.0 

7.0  MOS 

300  ±  96 

10 

80 

60 

20 

16/4/70 

39.15 

30.24 

5.4 

5.2  MOS 

4.4  ±  0.9 

7 

23/4/70 

39.45 

28.95 

5.2 

5.3  MOS 

3.8  ±  2.7 

6 

12/5/71  (06:25) 

37.68 

29.81 

5.5 

5.9  CGS 

26  ±  13 

13 

40 

20 

11 

12/5/71  (10:10) 

37.45 

30.13 

5.5 

3.5  ±  1 

14 

12/5/71  (12:57) 

37.61 

29.71 

5.4 

5.2  CGS 

10  ±  3 

14 

25/5/71 

39.20 

29.85 

5.8 

5.5  CGS 

12  ±  4 

13 

14/3/72 

38.79 

31.46 

5.4 

4.9  CGS 

3.4  ±  1.7 

8 

Brune's14  method  for  computing  seismic  slip  from  summed  moment  values  in  a  particular 
fault  zone  yields  a  value  of  1.3  cm/year  for  the  extension  in  Western  Turkey  over  1964-73.  This 
constitutes  reasonably  good  agreement  with  McKenzie's  estimate  of  2  to  3  cm/year  for  the  rel¬ 
ative  plate  motions.  Although  creep  appears  to  contribute  most  of  the  slip  in  other  parts  of  the 

1 3 

Mediterranean  and  Middle  East,  in  Western  Turkey  (at  least  over  the  last  ten  years)  it  seems 
to  have  been  taking  place  in  earthquakes. 


R.G.  North 
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Fig.  M-6.  Dependence  of  distortion  of  ray  paths  on  spreading  rate.  All  rav  oaths 
are  from  surface  events  on  ridge  axes.  All  rays  remain  in  vertical  plane  perpen- 
no  verf  i  g6  DiStE^ces  along  earth's  surface  are  marked  in  kilometers- 

0“to  ±70°al  eXaggeratlon-  Rays  are  spaced  by  5“  in  initial  takeoff  angle,  from 
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IV.  FOCAL  DEPTH 


A'  VIA  MAX, MUM 

The  shape  of  the  waveform  and  earl,  coda  contained  in  the  first  few  seconds  of  the  short- 
pertod  seismogram  ma,  be  determined  no,  onl,  b,  the  direr,  P-wa,e  bn,  also  bv  multipathed 
r-waves.  y  near-sonrce  and  race, ver  reverberation,  and,  the  case  of  shallow  events  by 
ne  uepth  phases.  The  contribution  of  each  of  these  secondary  arrivals  mav  be  nearly  the  same 
as  tha,  of  the  dtree,  wave.  Knowledge  of  the  onset  time  and  shape  of  the  later  arrivals  (and 
consequently,  their  removal,  is  often  necessary  to  obtain  an  accurate  representation  o,  the  ' 
sou  ce  function  „,e  short-period  band.  Determination  of  the  delay  time  of  pP  provides  the 
most  accurate  value  for  hypocentral  depth.  Assuming  that  these  secondary  arrivals  differ  in 
spectral  shape  from  the  direct  arrival  by  only  a  scale  factor  and  an  amount  due  to  differential 

— - ,he,r  de,ec,ion  may  *  acc°mpiishe<i  <=«>••«>  *„<,*  ^  echoes 
sample  of  thT'"  Wa!'el,'rm  °f  ‘he  Prlmar)'  a‘8nal-  «“  >°*  =Pectrum  contains  only  a  short  data 

Entropy  (M  F  I ZZT7  'h"'  'h'  a*  the  onsel  time  °r  the  echo.  Burg's  Maximum 

resZboL  error  spectra,  estimator  has  been  used  to  obtain  the  necessary 

The  value  of  using  the  M.  E.  method  for  computing  the  power  spectrum  of  truncated  real 
Sinusoids  has  been  demonstrated  by  Ulrych.'  cepstra,  analysis,  the  harmonics  are  complex 
and  o,  the  form  exp  [sqr,(- , )  .  w/w  |.  The  remarkable  resolution  obtained  by  using  Burg's 
,  chmque  ,s  shown  Fig.  ,V-t.  The  6-poin,  complex  M.  E.  prediction  error  operator  was  con- 
s  rue  ed  us.ng  cycles  of  the  complex.!, armonic  plus  50-percent  white  noise.  T  he  operator 
was  evaluated  around  the  uni,  circle  (at  10  times  the  sampling  interval).  The  resulting  spec- 
rum  p  otted  on  a  log  scale  shows  the  narrow,  well-defined  peak.  The  results  are  just  as  good 

ule  XT  r  r7“  °“e'third  °f  ‘  C,C'e  ““  Wi’h  -  mUCh  »  «•■»—  white  noise. 

ec  mque  in  cepstral  analysis  is  illustrated  in  Fig.lV-2.  The  input  seis- 

~ar;;:red  ? a  sourcc  fonc,ion  added  ,p ari"  -*■.  *  -  *  *  and 

,  a  a  "S  are  ,he  ™PUlSe  “f  the  LASA  short-period  instrument 

convolved  with  attenuated  impulses  with  travel  time  to  Q  ratios  If)  of  0.25  and  0.  The  log 

spec  rum  was  obtained  using  the  fas.  Fourier  transform.  The  phase  (imaginary  par,  of  the  log 

rum  was  unwrapped  and  a  linear  component  removed.  The  log  spectrum  of  the  instrument 
was  su  rac  e  in  the  band  (0.5  to  6  Hz)  where  sufficient  energy  exists  to  keep  the  phase  well- 
defined.  An  8- point  complex  M.  E.  prediction  error  operator  was  constructed  on  the  resting 
complex  frequency  series.  The  power  spectra  estimated  from  the  operators  show  well-defined 
peaks  detecting  the  echoes  tha,  are  no,  apparent  in  the  seismograms. 

Th,s  study  is  continuing  with  the  application  of  the  method  to  real  data. 

T.  E.  Landers 


B. 


FOCAL  DEPTH  FROM  SURFACE  WAVES 


recent  years,  considerable  attention  has  been  paid  to  the  possibility  of  using  surface- 
wa  e  spectra  to  determine  the  foca,  depth  o,  earthquakes.  This  method  has  obvious  application 
.scr,m, nation  the  case  of  shallow  events  for  which  depth  phases  are  unavailable  and  would 

f " couid  °niy  be  sh°-  -  -  ^ *■>*  -  *  -  -  excess  o,:r 
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Surface-wave  spectra  depend  in  a  complicated  manner  upon  the  source  mechanism,  focal 
depth,  and  elastic  parameters  of  the  earth.  Studies  to  determine  the  focal  depth  from  such 
information2,3  gener  ally  involve  some  form  of  fitting  scheme  to  observed  spectra  to  determine 
optimum  values  of  the  angles  describing  the  orientation  of  the  fault  plane  and  its  depth.  How¬ 
ever,  little  attention  has  been  paid  to  the  effect  of  imperfect  knowledge  of  earth  structure  on 
such  determinations,  and  generally  a  standard  continental  model  such  as  that  of  Gutenberg  and 
Bullen  has  been  used.  A  preliminary  study  has  been  made  of  the  dependence  of  spectral  ampli¬ 
tudes  of  Rayi'dgh  waves  on  reasonable  variations  in  shear-wave  velocity  structure.  The  ef  ect 

of  va.iable  eompressional  wave  velocities  and  densities  is  less  marked. 

Initial  studies4  to  determine  focal  depths  from  surface  waves  concentrated  on  the  period 
range  10  to  50  sec,  but  subsequent  work2  has  shown  that  waves  of  periods  less  than  ~  20  sec 
may  be  unprcdictably  affected  by  near-surface  variations  in  Q  and  should  not  be  used.  The 
amplitudes  of  shorter-period  waves  decrease  with  increasing  focal  depth,  but  it  is  difficult  to 
separate  this  effect  from  that  of  attenuation.  For  many  types  of  mechanism,  however,  the  theo¬ 
retical  amplitude  spectra  demonstrate  a  "hole"  at  periods  which  increase  with  increasing  source 
depth.  Unless  there  are  layers  of  low  Q  over  certain  depth  ranges  (which  is  certainly  the  case 
in  oceanic  structures  where  they  may  be  introduced  by  layers  of  sediment),  such  holes  will  not 
result  from  Q-variations.  The  effect  of  shear-wave  velocity  variations  on  the  periods  at  which 
holes  corresponding  to  particular  source  depths  occur  has  thus  been  chosen  as  suitable  for  stu  y. 

Crustal  structure  over  a  limited  area  can  often  be  fairly  well  determined  from  refraction/ 
reflection  experiments,  but  upper-mantle  structure  is  generally  inaccessible  by  these  methods. 
Normally,  upper-mantle  structure  (50-  to  500-km  depth)  is  obtained  from  surface-wave  disper¬ 
sion  curves,  whose  inversion  is  nonunique  and  provides  a  set  of  possible  models  satisfying  the 
data.  Some  constraints  on  these  models  occasionally  may  be  provided,  in  the  case  of  shear- 
wave  velocity  structure,  by  travel-time  data.  A  recent  study  of  upper-mantle  structure  un  er 
the  United  States  by  "hedgehog"  inversion  of  Rayleigh-wave  dispersion  curves  has  provided  sets 
of  models  satisfying  such  curves  determined  over  paths  of  400-  to  2000-km  length.  Rayleigh- 
wave  amplitude  spectra  for  sources  of  various  orientation  and  depth  within  these  models  ave 
been  computed,  and,  even  for  models  which  all  satisfy  the  same  dispersion  data,  the  variation 
of  the  position  of  the  spectral  hole  with  period  for  a  particular  source  depth  is  considerable. 

The  models  studied  to  date  are  those  determined  by  Biswas  and  Knopoff  to  satisfy  their  dis¬ 
persion  curves  R-3  (Western  U.  S.),  R-6  (North-Central  U.  S.),  and  R-ll  (South-Central  IU  ..). 

The  models  vary  shear-wave  velocities  in  three  layers  only,  designated  a  lid,  a  low-velocity 

zone,  and  a  layer  beneath  the  IV Z.  ,  ,  .. 

In  the  case  of  path  R-3.  f  models  determined  include  several  with  a  thin  high-velocity 

(4.6  to  4.95  km/sec)  layer  <  -low  the  Moho.  Amplitude  spectra  for  some  of  these  models 
show  spectral  holes  who.  pmu  :  arc  virtually  invariant  with  depth  over  i  5  to  50  km,  but  home 
of  these  structures  would  almost  certainly  be  excluded  by  shear-wave  travel-time  data.  Monels 
satisfying  paths  R-ll  and  R-6  are  more  reasonable,  but  even  for  these  the  variation  with  dept 
of  the  period  of  the  spectral  hole  differs  considerably  from  model  to  model.  For  all  models  for 
R-ll  there  is  negligible  difference  between  spectra  for  sources  at  40-  and  50-km  depth,  an 
the  amplitude  difference  could  easily  be  attributed  to  a  difference  in  seismic  moment.  For  some 
of  the  R-6  models,  there  is  a  similar  lack  of  distinction  between  spectra  for  depths  of  30,  40, 
and  50  km.  Spectra  for  a  vertical  strike-slip  fault,  observed  at  an  azimuth  of  30»  with  respect 
to  the  fault  strike,  for  depths  of  15,  20,  25,  30,  40,  and  50  km  for  Model  25  (R-6)  of  Biswas 


WB*wwwewwf 


and  Knopoff5  are  shown  in  Fig.  IV-3.  This  mode!  appears  he  in  no  way  extreme,  and  for  such 

a  fault  the  position  o,  the  hoies  does  no,  change  with  azimuth;  thus,  more  observations  would 
not  improve  the  depth  resolution. 

In  Table  IV-1.  the  period  at  which  the  spectral  hole  occurs  is  given  for  a  vertical  strike- 
P  fault  at  various  depths  in  seven  models  satisfying  dispersion  curve  R-6.  If  one  bears  in 
min  that  a  spectral  hole  cannot  be  observationally  resolved  to  better  than  a  *2-sec  period,  it 
can  e  seen  that  a  hole  at,  for  example,  36  sec  yields  a  depth  of  25  to  50  km.  Figure  IV- 4  shows 
the  upper-mantle  shear-wave  velocity  structure  for  these  seven  models. 


TABLE  IV-1 


PERIODS  (SECONDS)  AT  WHICH  SPECTRAI  wnr  v 


Depth 

(km) 

Model 

226 

23 

79 

179 

25 

213 

99 

Periods  (sec) 

15 

20 

22 

22 

21 

22 

23 

21 

20 

28 

28 

28 

27 

28 

28 

28 

25 

34 

34 

33 

32 

32 

32 

33 

30 

36 

36 

36 

36 

36 

36 

35 

40 

40 

38 

40 

38 

36 

36 

40 

50 

43 

38 

40 

42 

38 

36 

40 

Since  for  each  set  of  models  satisfying  a  particular  dispersion  curve  Biswas  and  Knopoff 
have  vaced  only  the  upper- man,, e  structure,  and  maintained  the  crustal  structure  cons  1,  „ 

,s  no,  surprising  tha,  there  is  little  variation  the  period  of  the  hole  for  depths  o,  5  to  zTkm 
owever,  it  is  certain  that  crustal  structure  is  subject  to  even  more  variation  than  that  of  the 
upper  mantle,  and  thus  the  resolution  of  me  method  for  events  shallower  than  20  km  is  probably 
even  poorer  than  demonstrated  here  for  those  in  excess  of  this  depth.  „  is  proposed  to  extend 
this  study  to  include  the  effect  of  variations  in  crustal  shear-wave  velocities 

in  theW?„rrofS!rn  ^  “  area  Wher'  8°°d  d8tB  °n  <h'  -.liable 

the  form  of  dispersion  curves,  the  accuracy  of  the  surface-wave  method  for  determination  of 

c  depth  ,s  limited.  In  regions  such  as  Central  Asia  where  the  structure  is  very  poorly  known 

the  resolution  will  deteriorate  further,  without  even  taking  into  account  the  probability  that  depth 

may  bo  introduced  spur,ousiy  by  pr°~  —  **  — -  r 

R.  G.  North 


28 


WRPWPPP 


1  I"  l|ll  MIlKlWIUgL..  . 


mirmm 


jw.uu  .419  wmtuMumg 


| 

E  i 


5, 

I 


REFERENCES 


1.  T.  J.  Ulrycli,  "Maximum  Entropy  Power  Spectrum  of  Truncated 
Sinusoids,"  J.  Geophys.  Res.  77,  1396  (1972). 

2.  Y-B.  Tsai  and  W-W.  Shen,  "Utility  of  Tsai's  Method  for  Seismic 
Discrimination,"  Texas  Instruments  Report  (30  November  1972). 

3.  D.  J.  Weidner  and  K.  Aki,  "Focal  Depth  and  Mechanism  of  Mid- Ocean 
Ridge  Earthquakes,"  J.  Geophys.  Res.  78,  1818  (1973). 

Y-B.  Tsai  and  K.  Aki,  "Precise  Focal  Depth  Determination  from 
Amplitude  Spectra  of  Surface  Waves,"  J.  Geophys.  Res.  76,  5729  (1970). 

N.  N.  Biswas  and  L.  Knopoff,  "The  Structure  of  the  Upper  Mantle 
Under  the  United  States  from  Dispersion  of  Rayleigh  Waves,"  Geophys. 

J.  R.  Astr.  Soc.  36,  515  (1974). 


4. 


5. 


—TO.2  r  — 


M  E.  PREDICTION 
ERROR  OPERATOR 

-  REAL  PART 

- IMAGINARY  PART 


Fig.  IV- 1.  Power  spectral  estimates  computed  from  6-point  maximum 
entropy  prediction  error  operator  designed  on  2  cycles  of  complex 
harmonic  with  50-percent  white  noise. 
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V.  SURFACE-WAVE  STUDIES 

A,  PREDICTING  RAYLEIGH -WAVE  AMPLITUDES  FROM  REFRACTION  MODELS 

Small-scale  lateral  variations  in  surface-wave  phase  velocity  cause  amplitude  and  directional 
anomalies  complicating  the  determination  of  M  which  is  critical  in  discrimination.  In  a  pre- 
vious  SATS,  Julian  reported  on  a  phase  velocity  model  for  the  western  U.S.  which  had  some  suc¬ 
cess  in  predicting  angles  of  incidence.  This  model  has  now  been  extended  into  the  Pacific  Ocean 
and  improved  through  the  inclusion  of  more  crustal  provinces.  Rays  traced  through  this  model 
from  LASA  or  an  explosion  site  show  patterns  of  focusing  and  defocusing  which  can  be  used  to 
predict  amplitudes. 

Models  have  been  developed  for  many  periods  between  15  and  40  sec.  The  improved  phase 
velocity  model  (Fig.  V-l)  shows  that  15-sec  structures  along  the  Pacific  Coast  are  much  more 
interesting  than  those  well  within  the  continent.  Generally,  the  ocean  is  about  0.16  km/sec  faster 
at  this  period,  with  the  transition  spread  over  some  200  km  and  distorted  by  local  anomalies. 

The  major  high  in  the  Oregon  area  is  caused  by  thin  crust  under  the  Columbia  Plateau  and  con¬ 
firmed  by  a  study  of  angles  of  incidence  of  Rayleigh  waves  at  LASA. 

When  rays  are  traced  through  this  model  and  arrive  at  LASA,  refraction  angles  of  30 “are 
common,  with  a  maximum  of  75°  Energy  should  be  focused  in  coming  across  the  northern  Cal¬ 
ifornia  coast  and  uefocused  in  the  Oregon-Washington  area.  Also,  it  is  apparent  that  multipath- 
ing  may  begin  at  the  coast,  if  not  before. 

As  a  test  of  this  model,  Rayleigh  amplitudes  at  various  stations  for  the  Faultless  explosion 
in  Nevada  were  predicted  by  ray  tracing.  The  separation  of  adjacent  rays  is  taken  to  be  inversely 
proportional  to  the  energy  concentration  in  the  wavefront.  Corresponding  data  were  obtained  by 
reading  the  maximum  vertical  long-period  displacements  from  records  of  WWSSN  and  LRSM  in¬ 
struments.  These  data  must  be  corrected  for  geometrical  spreading,  instrument  response,  atten¬ 
uation,  possible  double  couple  source  components,  and  the  theoretical  source  spectrum.  The 
results  (Fig.  V-2(a-c)]  show  an  accurate  prediction  of  the  high  amplitudes  in  the  Pacific  Northwest 
and  the  adjacent  lows  at  Berkeley  and  Bozeman.  Generally,  the  observed  amplitude  variations 
are  larger  than  those  predicted,  and  this  is  probably  because  no  attempt  has  yet  been  made  to 
correct  for  path -dependent  dispersion  effects. 

Thus,  we  find  that  lateral  refractions  of  over  30°  must  be  expected  for  short -period  Rayleigh 
waves  in  the  western  U.S.  These  should  mainly  be  caused  by  the  structure  near  the  Pacific  Coast. 
Ray  tracing  with  models  can  predict  necessary  Mg  corrections  of  up  to  0.5  unit  which  may  be  im¬ 
portant  in  seismic  discrimination. 

P.  Bird 
B.  R.  Julian 

B.  RAYLEIGH -WAVE  DISPERSION  FOR  THE  INDIAN  OCEAN 

To  gain  further  understanding  of  the  structural  evolution  of  the  oceanic  lithosphere,  a  project 
involving  the  study  of  regionalized  Rayleigh-wave  dispersion  for  the  Indian  Ocean  is  in  progress. 
Ihe  approach  of  our  regionalization  is  similar  to  that  carried  out  by  Forsyth‘S  on  the  Nasca  Plate. 

In  his  work,  Forsyth  used  magnetic  and  bathymetic  data  in  order  to  partition  the  spreading  lith¬ 
osphere  into  lateral  age  zones.  For  example,  if  we  model  the  lithosphere  with  three  oceanic  age 
zones,  a  suite  of  dispersion  curves  for  Rayleigh  waves  traveling  through  these  zones  may  be  used 
to  determine  phase  velocities  within  each  zone  by  requiring  that  the  sum  of  the  square  of  the  time 
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residuals  (expected  arrival  time  minus  observed  arrival  time)  be  a  minimum.  Furthermore, 
we  may  test  using  statistics  whether  a  particular  regionalization  model  with  three  age  zones  is 
significantly  better  than  a  model  with  just  two  age  zones.  For  the  study  on  the  Nasca  Plate, 
borsyth's  final  regionalized  model  required  three  oceanic  age  zones  and  two  continental  zones 
and  incorporated  anisotropy.  It  is  believed  that  by  studying  the  Indian  Ocean  in  a  similar  manner, 
some  important  questions  about  the  evolution  of  very  old,  oceanic  lithosphere  may  be  answered. 

To  date,  our  work  has  concentrated  on  building  up  a  suitable  data  base  of  dispersion  curves. 
Phase  velocities  for  vertical  component  Rayleigh  waves  originating  on  the  Mid -Indian  Ocean 
Ridge  system  have  been  determined  for  periods  16  to  150  sec  using  the  single -station  technique. 
This  technique  requires  that  the  initial  phase  at  the  source  be  known.  If  the  source  mechanism 
and  depth  for  an  earthquake  are  known,  we  may  calculate  the  initial  phase  using  the  theory  of 
Saito.  Thirteen  events  have  been  studied,  and  their  source  mechanisms  and  depths  have  been 
found.  From  these  13  events,  29  dispersion  curves  have  been  computed  for  the  paths  shown  in 
Fig.  V -3.  Total  random  errors,  such  as  those  associated  with  the  uncertainty  in  origin  time, 
digitization,  fault  plane  solution,  etc.  result  in  about  1 -percent  error  in  the  phase  velocity 
measurements. 

We  may  report  various  trends  in  these  data.  In  general,  wavetrains  from  the  Indi»-  ,^ean 
are  complicated  as  a  result  of  various  interference  effects  which  may  be  due  i-  to  strong 
lateral  variations  in  structure  which  are  known  to  exist  in  this  ocean  -»ong  attenuation  of  seis¬ 
mic  waves  is  observed  and  attributed  to  thick  sedimentary  la;  We  have  tentatively  correlated 

some  variations  of  thes*>  "path  averaged"  data  to  the  age  of  the  ocean  floor. 

H.  Patton 
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Fig.  V-3.  Rayleigh-wave  dispersion  paths  in  the  Indian  Ocean. 
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VI.  EARTH  HETEROGENEITY 


A.  AMPLITUDE -DISTANCE  CALIBRATION  CURVE 
FROM  DEEP  FOCUS  EARTHQUAKES 

The  reliability  of  seismic  magnitude  determinations  is  limited  ultimately  by  our  ability  to 
correct  for  the  effects  of  propagation  through  the  earth.  Magnitude  is  conventionally  determined 
from  a  relationship  of  the  form 

mb  =  log  (A/T)  +  B(h,  A)  (Vl-1) 


where  A  is  the  ground  amplitude,  T  is  the  corresponding  period,  h  is  the  focal  depth,  and  A  is 
the  epicentral  distance.  B(h,  A)  is  the  correction  factor  for  geometrical  spreading  and  attenu¬ 
ation  in  a  spherically  symmetric  earth.  The  values  of  this  function  given  by  Gutenberg  and 
Richter  are  still  the  accepted  standard,  although  refinements  have  been  proposed  by  Nuttli2  and 
by  Booth,  Marshall  and  Young,  among  others.  We  propose  here  another  refinement  based  on 
some  500  P-amplitude  data  for  14  deep  events  in  different  Benioff  zones  around  the  world.  The 
impulsive  nature  of  P-signals  from  deep  focus  earthquakes  as  well  as  the  avoidance  of  near¬ 
source  upper-mantle  heterogeneity  make  them  well-suited  to  studying  propagation  through  the 
lower  mantle. 


The  amplitude  data  were  corrected  for  instrument  response  and  the  effects  of  reflection/ 
transmission  at  the  free  surface,  the  Moho,  and  the  650-km  discontinuity.  Correction  for  the 
radiation  pattern  was  also  made  using  fault  plane  solutions  determined  from  P  first-motion  data 
and  S-wave  polarization  angles.  All  amplitude  data  requiring  a  radiation  correction  in  excess  of 
a  factor  of  3.0  were  deleted  from  the  study.  A  correction  for  the  source  spectrum  was  made 
assuming  that  the  amplitude  spectrum  has  a  corner  frequency,  below  which  it  is  flat  and  above 
which  it  decreases  according  to  a  power  law.  Following  Wyss  and  Molnar,4  we  assumed  that  for 


the  events  in  this  study,  the  corner  frequencies  are  above  0.2  Hz,  and  included  only  those  am¬ 
plitude  data  for  which  the  measured  periods  were  greater  than  5  sec.  The  averages  of  these 
amplitudes  for  all  events  were  then  used  to  equalize  the  moments  to  a  common  value.  Through 
a  process  of  iteration,  we  then  determined  the  average  values  of  amplitude  in  2°  cells  and  the 
station  anomalies  at  different  stations. 

Figure  Vl-1  shows  our  amplitude  data,  averaged  in  2°  cells,  after  correction  for  station 
anomalies.  These  amplitude  data  roughly  correspond  to  a  source  depth  of  550  km.  The  solid 
curve  if  the  ray  theoretical  amplitude  from  a  model  which  satisfies  the  P  travel-time  data  of 
these  deep  earthquakes.  The  sawtooth  appearance  of  ray  theoretical  amplitude  is  merely  a  char¬ 
acteristic  of  the  amplitvdc  calculation  from  ray  theory.  The  dashed  curve  is  our  version  of  the 
amplitude  distance  curve  after  smoothing  the  ray  theoretical  amplitudes  to  correspond  to  the 
data  as  closely  as  possible. 

Figure  Vl-2  shows  our  amplitude-distance  curve  in  comparison  with  the  Gutenberg  and 

j 

Richter  curve.  The  shapes  of  the  two  curves  differ  considerably,  possibly  due  to  frequency- 
dependent  attenuation,  as  the  Gutenberg-Richter  curve  was  determined  from  short-period  P- 
amplitudes  whereas  our  data  have  periods  greater  than  5  sec.  Incidentally,  our  curves  have 
more  similarity  with  those  determined  by  Booth,  Marshall  and  Young3  than  with  the  one  deter- 
mined  by  Nuttli.  Currently,  the  effect  of  this  curve  upon  the  reliability  of  m.  determinations 


is  being  evaluated. 


M.  Sengupta 
B.  R.  Julian 


37 


B. 


MORE  UP'.-  ER-MANTLE  REFLECTIONS 


In  the  last  SATS.5  we  reported  on  the  observation,  at  LASA,  of  intei  mediate,  short-period 
phases  between  P  and  pP  from  deep  earthquakes  in  western  South  America.  An  interpretation 
that  may  be  given  to  these  phases,  which  are  called  pdP,  is  that  they  are  reflections  from  bound¬ 
aries  within  the  upper  mantle  above  these  sources;  alternate  interpretations  are  that  they  are 
aftershocks  or  reflections  from  contrasts  beneath  the  array.  In  order  to  test  these  various 
hypotheses,  we  have  extended  the  study  to  other  zones  where  deep  seismicity  exists.  The  results 

°‘  thiS  eXtenSi°n  are  ShOWn  in  FiS-  VI~3  the  time  difference  pP-pdP  is  plotted  as  a  function 

of  distance  for  the  five  regions  of  deep  seismicity  studied.  We  point  out  that  a  pdP  phase  is  not 
clear  from  every  deep  earthquake  studied  at  LASA.  and  only  the  data  from  unambiguous  arrivals 
are  plotted  in  Fig.  Vl-3.  The  solid  lines  running  across  that  figure  represent  the  depth  of 
source  from  standard  pP-P  tables.  Although  no  completely  clear  pattern  is  developed  most 
of  the  arrivals  seem  to  indicate  a  depth  of  the  reflection  point  between  200  and  300  km  Since 
the  upper-mantle  structure  above  these  deep  sources  is  complicated  by  cool  lithospheric  slabs 
descending  into  the  mantle,  these  reflections  may  not  be  given  to  any  simple  interpretation  in 
terms  of  plane,  horizontal  boundaries.  There  is  one  rather  loose  correlation  that  may  be  rel¬ 
evant.  The  seismicity  of  some  of  the  deep  earthquake  zones  is  characterized  by  a  dearth  or 
gap  in  the  activity  at  intermediate  depths,  generally  between  300  and  500  km.  This  gap  is  partic¬ 
ularly  well-defined  beneath  South  America  and  the  New  Hebrides,  less  so  beneath  the  Sea  of 
Okhotsk,  and  only  vaguely  apparent  beneath  the  Bonin  and  Tonga  Islands.  In  Fig.  Vl-4,  we  have 
plotted  the  implied  reflection  depths  of  these  pP  precursors  along  with  the  extent  of  the  seismic¬ 
ity  gaps,  if  clear,  in  the  five  regions  studied.  From  this  figure  it  appears  that  many  of  the  re¬ 
flection  points  fall  near  the  top  boundary  of  the  intermediate  zones  of  low  or  no  seismicity  If 
this  correlation  is  a  physical  one.  it  would  imply  rather  sharp  changes  in  the  elastic  parameters 
of  the  upper  mantle  associated  with  the  upper  boundaries  of  the  aseismic  zones.  Since  the  pP 
precursors  seem  to  be  fairly  common,  we  can  reject  the  aftershock  hypothesis  with  some  confi¬ 
dence.  We  likewise  tend  to  reject  a  source  beneath  the  array  as  a  cause  of  these  phases  since 
they  are  not  found  following  the  short-period  P-waves  of  shallow  events,  however  the  clarity  of 
such  a  reflection  may  depend  on  the  frequency  content  of  the  primary  wave. 

We  would  like  to  correct  two  errors  in  our  previous  discussion5  of  these  phases.  The  array 

traces  in  bigs.  lll-8(a)  and  (b)  of  Ref.  5  are  mislabeled.  They  should  read,  from  the  top  Bi  F3 
F4,  A0,  B3.  C4,  B4.  Cl.  C2.  B2,  C3,  and  D3.  It  has  been  pointed  out  to  us6  that  the  phrase 
"plane  of  constant  latitude,"  used  in  the  text  and  in  the  caption  to  Fig.  111-10,  is  not  correct.  A 
surface  of  constant  latitude  is  a  cone,  not  a  plane.  The  seismicity  discussed  was  in  fact  pro¬ 
jected  onto  the  plane  of  a  great  circle  which  dips  25“  with  respect  to  the  equatorial  plane  at  a 
direction  of  70  °W  longitude  on  that  plane 

M.  Lin 

J.  R.  Filson 


C.  PRECURSORS  TO  S  AND  SKS  FROM  DEEP  FOCUS  EARTHQUAKES 

LASA  data  have  been  searched  for  arrivals  preceding  the  S-  and  SKS-phases  which  are  ex¬ 
pected  to  result  from  the  conversion  of  shear  wave  energy  to  compressional  wave  energy  at 
upper-mantle  discontinuities.  Numerous  authors  have  used  precursors  to  PP  and  P'P’  in  study¬ 
ing  such  discontinuities.  The  purpose  of  this  study  is  to  determine  the  local  structure  of  upper- 
mantle  discontinuities  under  seismic  array  sites,  information  about  which  could  be  valuable  for 
many  studies  involving  array  data. 
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The  criteria  used  in  identifying  a  precursor  to  a  main  phase  are  slowness,  azimuth,  and 
amplitude.  The  slowness  and  azimuth  of  an  arrival  are  found  by  picking  the  arrival  time  of  the 
phase  on  as  many  of  the  array  channels  as  possible  and  best  fitting  a  plane  wave  front.  A  pre¬ 
cursor  resulting  from  the  conversion  of  shear  to  compressional  energy  will  have  the  same  slow¬ 
ness  and  azimuth  as  the  main  shear  arrival.  The  Velocity  Spectral  Analysis  (VESPA)  process 
is  currently  being  applied  to  determine  the  slowness  of  arrivals  preceding  the  S-  and  SKS-phases. 

Events  recorded  at  LASA  having  a  magnitude  of  5.4  or  greater  and  a  depth  greater  than 
200  km  were  selected  for  this  study.  S-to-P  conversion  is  expected  only  beyond  critical  distances 
of  82“  for  the  420-km  discontinuity,  and  89“  for  the  670-km  discontinuity.  Events  were  selected 
for  this  study  with  distances  greater  than  89“  so  that  precursors  to  the  S-phase  converted  at  both 
discontinuities  could  potentially  be  seen  for  each  event.  The  shadow  zone  due  to  the  core  reduces 
the  amplitude  of  P-  and  S-phase  arrivals  for  events  with  distances  greater  than  102“  so  that  no 
events  with  distances  greater  than  102“  were  selected. 

The  S-  and  SKS-phases  were  chosen  for  the  study  because  of  their  low  values  of  slowness 
and  their  arrival  in  a  portion  of  the  record  where  there  are  comparatively  few  arrivals  of  other 
phases.  The  short-period  vertical  component  and  the  three  long-period  components  were  ex¬ 
amined  for  each  of  the  26  events  selected. 

Precursors  were  most  commonly  seen  in  the  range  60  to  70  sec  before  the  S-  and  SKS-phase 
arrivals.  The  slowness  values  of  the  S-phase  are  quite  different  from  those  of  the  SKS-phase, 
so  that  the  corresponding  main  phase  of  a  given  precursor  can  be  clearly  determined.  Precur¬ 
sors  were  most  prominent  on  the  vertical  components.  Precursors  to  the  S-  and  SKS-phase 
arrivals  generated  at  the  670-km  discontinuity  are  expected  about  65  sec  before  the  main  phase, 
and  for  the  420-km  discontinuity  about  45  sec  before  the  main  phase. 

The  precursors  varied  from  a  single  peak  to  a  coda  of  several  wavelengths.  The  period  of 
the  precursors  averaged  around  8  to  10  sec  on  the  long-period  data,  and  around  1  sec  on  the 
short-period  data. 

Figure  VI-5  shows  a  precursor  to  the  S-phase  arriving  70  sec  before  the  main  phase  on 
long-period  channels.  A  precursor  to  the  SKS-phase  is  shown  in  Fig.  VI-6  arriving  64  sec  before 
the  main  phase  on  short-period  channels. 

The  spectrum  of  the  precursor  coda  is  being  analyzed  in  an  effort  to  gain  insight  into  the 
detailed  structure  of  the  transition  zones. 

C.  T.  Bolt 
B.  R.  Julian 

D.  SCATTERING  EFFECTS  ON  PcP/P  AMPLITUDE  RATIOS 

A  strong  source  of  amplitude  scattering  for  short-peviod  data  is  the  crust  and  upper-mantle 
structure  under  the  receiver.  This  is  impossible  to  estimate  with  a  single  station;  however,  an 
array  of  closely  spaced  seismographs  gives  a  clear  picture  of  the  complexity  of  scattering  effects. 
Studies  by  K.  Aki  and  J.  Capon  show  that  a  considerable  amount  of  amplitude  scattering  at 

LASA  can  be  caused  by  a  medium  with  small  random  fluctuations  of  velocity  of  a  percent  or  two 
about  the  mean. 

This  has  disastrous  consequences  for  studies  based  on  amplitude  ratios  of  short-period 
phases  such  as  P  ana  PcP.  To  illustrate  this,  six  events  from  a  recent  study9  were  chosen 
which  had  good  signal-to-noise  ratios  at  all  LASA  subarrays  with  distances  of  44°  to  63“  along 
a  narrow  azimuthal  sector  302“  to  317“.  Four  events  were  clustered  in  the  region  from  Andreanof 
Island  to  Amchitka  Island,  and  two  events  were  located  near  Kamchatka.  Table  VI-1  contains  a 
list  of  the  events  used. 
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PcP/p  amplitude  ratios  were  then  computed  for  each  subarray.  In  Table  VI-1,  the  mean 
Pc  P/p  ratio  and  standard  error  across  LASA  for  each  event  are  tabulated.  Note  that  the  stand¬ 
ard  errors  decrease  with  increasing  epicentral  distance  except  for  event  2.  This  is  deceptive, 
however,  because  the  mean  PcP/P  values  also  decrease  with  distance.  In  the  last  column  of 
Table  VI-1,  the  normalized  standard  errors,  i.e.,  standard  error/mean,  seem  to  be  more  or 
less  comparable  at  all  distances.  These  statistics  show  that  PcP/p  ratios  on  subarray  beams 
can  easily  vary  by  ±50  percent  relative  to  the  mean  PcP/P  ratio  across  the  LASA. 

A  measure  of  similarity  of  two  amplitude  patterns  is  given  by  their  correlation  coefficient. 
The  log  amplitudes  of  the  P  and  PcP  phases  were  computed  at  all  21  subarrays  for  the  six 
events.  The  mean  of  each  log  amplitude  pattern  was  removed  yielding  a  pattern  of  positive  and 
negative  log  amplitudes  for  each  phase.  Fo.'  event  K,  these  can  be  represented  by  the  vectors 


and 

PcPK  =  {PcP.K) 

l 

where  i  =  1,  21  corresponding  to  the  21  subarrays  of  LASA. 

The  correlation  coefficient  of  A  and  B  is  defined  as 

C(A,B)  =  ~=S== 
v  (A  A)  (BB) 

Three  sets  of  correlation  coefficients  were  calculated,  namely  C(PK  PcPK)  C(PK  PK+11 
K  K  +1  •  #  * 

and  C(PcP  ,PcP  ).  The  first  is  the  correlation  coefficient  of  P  and  PcP  patterns  from  event 

K;  the  second  involves  the  P  patterns  from  two  adjacent  events  in  Table  VI-1;  and  the  third  co¬ 
efficient  is  for  PcP  patterns  from  two  adjacent  events.  These  calculations  are  shown  in  Fig.  VI-7. 

For  the  events  at  distances  less  than  50%  P  and  PcP  patterns  are  essentially  uncorrelated, 
whereas  at  about  60%  the  patterns  are  moderately  correlated.  On  the  other  hand,  PcP  patterns 
from  adjacent  events  are  well  correlated,  as  shown  by  the  dashed  lines  between  adjacent  event 
distances,  except  for  the  correlation  between  the  events  at  47°  and  58%  the  largest  distance  gap. 
The  P  patterns  are  slightly  less  correlated  than  the  PcP  patterns,  but  generally  are  better  than 
the  correlation  between  P  and  PcP  patterns  for  each  of  the  adjacent  events. 

The  implication  of  these  data  is  that  P  and  PcP  amplitude  patterns  at  LASA  can  be  duplicated 
by  closely  located  events,  the  PcP  pattern  being  more  consistent  than  the  P  pattern  with  varia¬ 
tions  in  epicentral  distance.  This  result  indicates  that  the  amplitude  Fluctuations  of  subarray 
sum  data  are  dominated  by  complex  transmission  and  diffraction  of  the  seismic  signals  by  the 
structure  under  LASA  rather  than  near  the  source  region  or  deep  in  the  mantle.  Most  likely  this 
is  the  case  under  many  seismograph  stations,  and  suggests  that  single-site  measurements  of 
PcP/P  ratios  are  not  reliable  for  determining  parameters  of  the  core  mantle  boundary. 

By  forming  beams  of  the  P  and  PcP  phases,  the  incoherent  amplitude  fluctuations  are 
averaged  out,  yielding  more  stable  estimate  of  PcP/P  ratios.  If  the  PcP/P  ratios  measured 
at  different  subarrays  are  assumed  to  be  indt;  endent  with  equal  means  and  variances,  then  one 
can  expect  that  the  standard  errors  of  PcP/P  ratios  obtained  from  the  beams  are  1/V21  (~0.22) 
times  the  standard  errors  of  the  single  subarray  measurements.  From  Table  VI-1  this  implies 
that  the  standard  errors  for  each  ratio  using  beamed  data  are  about  10  percent  of  the  mean  value, 
which  is  a  more  acceptable  standard  error. 

C.  W.  Frasier 


E.  SHORT-PERIOD  CODA  OF  A  LOCAL  EVENT  AT  LASA 

It  is  widely  accepted  that  the  seismic  coda  of  small  local  earthquakes  is  largely  due  to 
scattering  by  inhomogeneitics  in  the  earth's  crust,  as  suggested  by  Aki.  The  actual  mechanism 
of  this  scattering  is,  as  yet,  unsatisfactorily  explained  and  must  be  understood  to  accurately  de¬ 
fine  those  parts  of  a  seismogram  due  to  source  effects  and  those  due  to  path  effects. 

It  is  possible  to  discern  at  least  two  processes  contributing  to  codas.  The  first  is  strong 
forward  scattering  discussed  by  Chernov.**  This  process  has  the  effect  oi  summing  a  large  num¬ 
ber  of  multipathed  signals.  Essentially,  all  the  energy  in  such  coda  travels  along  the  ray  path 
connecting  the  source  to  the  receiver.  Additionally,  as  a  wave  passes  a  receiver  it  may  be  scat¬ 
tered  locally  by  smaller  scale  inhomogeneities  in  the  vicinity  of  the  receiver.  This  phenomenon 
of  multiple  scattering  is  tire  second  process  contributing  to  seismic  codas.  According  to  the 
analysis  of  Ilowe,*2  which  predicts  an  equipartition  of  energy  between  all  wavenumbers,  the 
energy  in  this  section  of  the  coda  arrives  from  all  directions.  A  superposition  of  these  two  proc¬ 
esses  may  provide  a  good  qualitative  description  of  P-,  S-,  and  surface-wave  coda. 

Figure  VI-8  shows  the  energy  in  the  frequency  band  1.0  to  2.0  Hz  of  a  short-period  record 
at  the  F2  subarray  of  LASA  for  a  strip-mi  ling  blast  113  km  away,  as  a  function  of  time.  The 
peak  energy  in  this  band  is  due  to  the  arrival  of  the  S-wave.  The  high-resolution  wavenumoer 
spectrum  for  this  section  of  the  record,  shown  in  the  figure,  indicates  that  virtually  all  the 
energy  is  arriving  from  the  direction  of  the  blast.  The  wavenumber  spectrum  for  the  later  S-coda 
indicates  that  energy  is  arriving  from  all  directions  with  shear-  and  surface -wave  velocities. 

Figure  VI-9  shows  the  local  energy  in  the  frequency  band  0.2  to  1.0  Hz.  In  this  band,  the 
peak  energy  occurs  at  the  arrival  of  the  Rayleigh  wave.  The  wavenumber  spectrum  shows  that 
virtually  all  the  energy  in  this  section  of  the  record  arrives  from  the  direction  of  the  source, 
with  a  velocity  appropriate  to  2-sec  Rayleigh  waves  under  LASA.  The  coda  of  these  arrivals 
produces  a  wavenumber  spectrum  that  again  shows  that  energy  arrives  from  all  directions,  with 
velocities  ranging  from  P-  to  surface-wave  velocities.  These  results  bear  out  the  salient  points 
of  the  model. 

Before  this  model  can  be  effectively  applied  to  local  events  and  teleseisms,  it  is  necessary 
to  determine  more  accurately  the  nature  of  the  P-,  S-,  and  surface-codas.  This  information 
can  be  used  to  determine  values  for  the  correlation  distances  and  turbidity  coefficients  which 
characterize  the  type  of  scattering  in  the  crust  under  LASA. 

J.  Scheimer 
T.  E.  Landers 
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Fig.  VI-1.  Amplitude-distance  curve.  Circles  are  average  amplitude  (microns) 
in  2°  cells;  vertical  lines  are  standard  errors  of  the  average;  solid  line  is  ray 
theoretical  amplitude  from  model  MKS1(P),  satisfying  travel -time  data  from  deep 
focus  earthquakes.  Note  ray  theoretical  amplitude  has  been  adjusted  vertically 
to  satisfy  amplitude  data.  Dashed  line  is  smooth  amplitude-distance  curve. 


Fig.  VI-2.  Comparison  of  amplitude-distance  curves.  Solid  line  is  amplitudes 
(microns)  of  Gutenberg  and  Richter  for  earthquake  of  magnitude  6.0  and  period 
of  1  sec.  Dashed  line  is  our  version  of  amplitude -distance  curve. 
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Fig.  VI-3.  Time  interval  pP-pdP  at  LASA  from  deep 
sources  in  various  regions.  Three  solid  lines  repre¬ 
sent  depths  of  reflection  points  in  kilometers. 


REGION 


Fig.  VI-4.  Correlation  between  implied  reflection  point  depths 
(points)  and  seismicity  gaps  (arrows)  in  various  regions. 
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Fig.  VI-5.  Long-period,  subarray  sums  from  LASA  showing  precursor 

£  !'&'[, 22,?' SbelT  »tase  (™m  earthquake  occurring 

S  9  1  at  Z2h  52m  15s,  530  km  deep  beneath  20.3°S,  177.8°W. 


ShortTperl°d*  subarray  sums  from  LASA  showing  precursor 
HrmKnnP24STe  arr*'vl^^^>^  sec  before  main  phase  from  earthquake  occur - 
179  5°W  nUary  ^  1  33m  °3-5S>  595  km  deep  beneath  21  . 9 ° S, 


Fig.  VI-7.  Correlation  coefficients  of  P  and  PcP  patterns 
of  same  event  at  LASA,  of  P  patterns  from  adjacent  events, 
and  PcP  patterns  from  adjacent  events. 


Fig.  VI-8.  High-resolution  spectra  for  S-arrival  and  S-coda 
in  frequency  band  i.O  to  2.0  Hz  for  strip-mining  blast  near 


Fig.  vi-g.  High-resolution  spectra  for  Rayleigh-wave  arrival 
and  Rayleigh-wave  coda  in  frequency  band  0.2  to  1.0  Hz  for  same 
blast  in  Fig.  VI-8. 
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VII.  GENERAL  SEISMOLOGY 

A.  SEISMICITY  OF  CENTRAL  ASIA  AND  TECTONICS 

OF  THE  BAIKAL  RIFT  ZONE 

1  2 

In  two  earlier  reports  '  we  presented  fault-plane  solutions,  horizontal  slip-ve  :tors,  and 
paleopoles  for  Central  Asia.  We  also  presented  a  hypothetical  plate-model  for  Central  Asia. 

Figure  VII-1  shows  xir  seismicity  map  for  Central  Asia.  The  seismicity  is  plotted  for  the 
region  bounded  by  18  °N  and  58  °N  and  by  60  °E  and  140  “E.  Historical  as  well  as  instrumental 
seismicity  are  plotted  on  the  map.  The  historical  seismicity  is  obtained  from  the  ARPA  seis¬ 
micity  map.  The  instrumental  data  until  1961  are  also  obtained  from  the  ARPA  map.  Instru¬ 
mental  data  from  1962  to  1971  were  taken  from  the  PDE  cards.  All  magnitudes  are  >4.0  (m^h 
The  epicenters  are  shallow  (< 7 0  km)  except  a  few  in  the  Burma  region  and  those  of  the  Hindukush 
Zone  which  are  of  intermediate  depth  (between  70  and  300  km).  The  very  deep  earthquakes 
(~500  km)  near  Korea  are  interpreted  as  being  due  to  the  Pacific  plate  dipping  down  under  Asia. 

Even  though  the  seismicity  in  Central  Asia  is  diffused  over  a  wide  region,  we  are  still  able 
to  identify  zones  of  concentrated  seismicity.  It  is  along  these  zones  that  the  plate  boundaries 
were  drawn  in  Fig.  IV -12  of  Ref.  2. 

The  Baikal  Rift  Zone:-  This  zone  is  situated  on  the  boundary  of  the  Siberian  Platform 
and  exhibits  singular  features.  We  were  not  able  to  get  enough  data  to  do  more  than  two  focal 
mechanism  solutions  in  the  region,  since  most  of  the  earthquakes  between  1963  to  1972  were 
<5.5.  However,  Misharina  has  done  fault-plane  solutions  for  these  events  using  data  from 
FSSR  stations,  and  has  plotted  the  axes  of  tension  and  compression.  Some  of  the  faults  in  the 
region  are  shown  in  Fig.  15  of  Ref.  3.  All  earthquakes  in  the  lake  are  pure  normal.  Going  from 
the  northeastern  end  of  the  lake  toward  the  Stanovoy  mountains,  the  focal  mechanisms  are  still 
mostly  normal  but  some  have  a  large  component  of  strike  slip.  Our  solutions  Nos.  10  and  19 
(see  Ref.  1)  are  in  this  region.  No.  10  is  pure  normal,  wnile  No.  19  is  strike  slip  with  a  small 
component  of  thrust,  which  agrees  with  Misharina' s  solutions.  But  if  we  move  along  the  main 
Sayan  fault,  which  extends  at  an  acute  angle  of  40°  to  the  western  flank  of  the  rift  zone,  and 
which  has  no  depressions  of  the  rift  type,  all  the  ocal  mechanisms  are  of  purely  compressional 
type.  If  we  move  even  further  along  the  Si  berm. -W.  China  plate  boundary  (which  is  discussed 
later),  the  solutions  are  still  purely  compressional  -  Solutions  6,  18,  21,  and  23.  Misharina3 
indicates  the  region  where  the  sudden  change  from  extension  to  compression  takes  place.  But 
below  this  boundary,  the  earthquakes  within  the  plate  are  of  pure-strike  slip  type,  with  a  small 
component  of  thrusting  in  some  cases.  Solution  14  and  the  Gobi-Altai  earthquake^  solutions  are 
of  this  type.  The  level  of  seismicity  in  this  region  is  higher  than  in  the  rift  zone. 

Thus,  there  is  a  region  of  near-horizontal  stretching  at  the  upper  end  of  the  lake  and  near- 
horizontal  compression  beyond  the  lower  end  of  the  lake.  Florensov5  claims  that  there  are  signs 
of  increase  of  rift  depression  in  width  and  depth  and  of  active  rift  structures  in  the  Stanovoy 
mountains,  and  that  these  indicate  the  continuity  of  development  of  the  rift  zone  in  the  eastern 
direction.  The  rift  structures,  according  to  Florensov,  are  along  the  ancient  faults  of  the  Pre- 
Cambrian  and  early  Paleozoic  eras. 

The  rift  zone  is  characterized  by  an  intense  isostatic  instability  with  large  negative  gravity 
anomalies  which  some  claim  may  be  as  large  as  -100  mgals'.!  Borisov6  is  of  the  opinion  that 
beneath  the  entire  rift  is  a  very  deep  graben-like  structure  that  goes  down  into  the  upper  mantle 
and  is  filled  with  crust  mantle  mixture. 
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Lubimova  has  conducted  a  study  of  heat -flow  data  in  the  lake  and  on  the  land  surrounding 
the  lake.  It  is  obvious  that  the  Baikal  Zone  is  a  region  of  abnormally  high  heat  flow  (the  highest 
value  being  3.4  11FU).  But  the  region  just  beyond  the  southwest  tip  of  the  lake  shows  lower  values. 
Lubimova  and  Polyak  note  this  "clear  boundary"  separating  the  Baikal  Rift  Zone  from  the  "rel¬ 
atively  uniform  field  of  somew'hat  low'er  values  in  the  adjacent  part  of  the  Siberian  platform,  that 
of  the  Irkutsk  plateau  (~1.0  HFUL"  Note  that  the  transition  from  high  to  low  heat-flow  values 
occurs  at  the  exact  same  spot  as  the  transition  from  extensional  to  compressional  stresses. 

Thus,  we  find  a  most  unusual  phenomenon  occurring  on  the  boundary  between  the  Siberian 
and  the  W.  China  plates.  There  seems  to  be  compression  on  the  part  from  the  Hindukush  to  the 
southern  tip  of  Lake  Baikal,  but  from  there  to  the  Stanovoy  region  there  is  spreading.  The 
W.  China  plate  thus  seems  to  be  turning  in  a  clockwise  direction  relative  to  the  Siberian  plate 
which  is  assumed  fixed,  about  a  point  at  the  southern  tip  of  Lake  Baikal. 

S.  Das 

B.  PROPAGATION  OF  TI1E  NORSAR  SUBARRAY  TIME  CORRECTIONS 

INTO  ASEISMIC  REGIONS 

The  analysis  of  the  NORSAR  time  anomalies  has  reached  its  conclusion  with  the  implementa¬ 
tion  of  the  subarray  time  corrections  into  the  Analysis  Console  System  and  into  various  computer 
programs.  Before  the  corrections  were  implemented,  however,  it  was  necessary  to  extrapolate 
the  time  errors  into  regions  that  were  aseismic  or  lacked  sufficient  data  to  yield  a  value  for  the 
s  tation  correction.  The  process  used  for  the  extrapolation  of  the  data  is  the  same  one  success - 
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fully  used  earlier  on  the  LAS  A  station  corrections  for  core  phases. 

The  first  step  in  the  process  is  to  plot  the  subarray  time  errors  on  a  graph  of  azimuth  vs 
error  and  a  plot  of  distances  vs  error.  From  these  two  plots,  it  can  be  determined  what  type 
of  function  will  best  approximate  the  suharray  error.  In  the  case  of  the  NORSAR  subarrays,  the 
function  seemed  to  have  a  linear  term  or  a  quadratic  term  with  distance  and  was  definitely  a 
three-cycle  Fourier  series  type  function  with  azimuth.  A  good  example  of  this  is  subarray  03C 
and  is  shown  in  Fig.  Vll-2.  This  generally  represents  the  functional  relationship  of  the  NORSAR 
suharray  time  anomalies  with  distance  and  azimuth.  The  next  step  in  the  procedure  was  to  fit 
the  desired  function,  in  a  least-squares  sense,  to  each  set  of  subarray  time  residuals.  The 
function  that  was  used  had  21  coefficients,  and  the  fit  resulted  in  a  root-mean-square  error  that 
was  generally  about  0.07  sec.  The  evaluation  of  this  function  at  regions  that  are  aseismic  will 
yield  a  fair  value  for  the  subarray  station  correction  in  those  areas. 

However,  additional  refinement  of  the  functional  values  was  necessary.  Since  the  final  cor¬ 
rections  were  to  be  stored  in  the  form  of  a  table,  the  function  was  evaluated  at  each  cell  in  the 
table  (cell  size  is  10°  in  azimuth  by  0.2  sec/deg  in  dT/dA).  Where  actual  data  existed  at  a  cell 
point,  it  was  used  rather  than  the  function  value;  where  there  was  no  actual  data,  the  function 
value  was  used.  The  resulting  table  was  then  smoothed  somewhat  to  reduce  the  sharp  boundaries 
between  some  cells.  This  final  table  was  then  used  as  the  NORSAR  subarray  correction.  A 
graphic  example  of  the  final  output  for  subarray  03C  is  shown  on  the  polar  plot  in  Fig.  VII-3. 

Here,  the  radial  component  is  dT/dA  starting  at  4.0  in  the  center  and  extending  to  10.0  at  the 
outer  circumference. 

R.  M.  Sheppard 
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C.  SINGLE -CHANNEL  EVENT  DETECTOR  IN  REAL  TIME 

In  the  previous  SATS,^  a  single-channel  event  detector  was  described  which  operates  in  real 
time  on  a  PD  P-7  computer.  Here,  we  present  results  of  using  this  program  on  a  set  of  480  syn¬ 
thetic  events  derived  from  16  real  events  recorded  at  NORSAR. 

Briefly,  the  detection  scheme  compares  the  spectral  power  of  the  seismogram,  averaged 
over  a  short  time  window,  with  the  spectral  power  averaged  over  a  long  window  preceding  the 
short  window.  A  detection  occurs  when  the  short-term  power  (averaged  over  several  frequencies 
of  interest)  exceeds  the  long-term  power  (averaged  over  the  same  frequencies)  by  a  factor  of  T, 
the  threshold  parameter.  When  a  detection  is  declared,  the  long-term  average  (LTA)  is  fixed, 
while  the  short-term  average  (STA)  is  continuously  updated  in  time.  As  the  event  coda  decays 
in  time,  the  ratio  of  STA  to  LTA  decreases  until  it  falls  below  T,  at  which  time  the  detector 
turns  off.  If  no  new  detections  occur  within  an  arbitrary  time  lag  L,  then  the  LTA  is  updated, 
and  the  detection  scheme  continues. 

Experiments  have  shown'  that  the  frequencies  1.0,  1.5,  2.0,  and  2.5  Hz  were  suitable  for 
calculating  the  STA  and  LTA  for  NORSAR  data.  Other  parameters  used  were  a  threshold  T  =  2.5, 
a  lag  L  =  10  sec,  a  short-term  averaging  window  of  3.0  sec,  and  a  long-term  averaging  window 
of  30  sec. 

NORSAR  beams  of  16  events  were  formed  using  the  center  sensors  of  the  NORSAR  subarrayr. 
The  events  sample  different  azimuths  with  distances  from  38. 3 “to  77.8“  from  NORSAR.  Seven 
minutes  of  noise  from  5  NORSAR  channels  were  selected  from  13:02:00  to  13:19:00  GMT  on 
23  February  1972.  The  noise  channels  are  1A0,  1B0,  2B0,  3B0,  and  4B0.  Each  event  beam 
was  added  to  a  set  of  the  5  channels  at  1-min.  intervals  from  13:03  to  13:08,  scaled  down  by  a 
factor  of  0.5  or  drop  in  mfa  of  0.3  relative  to  the  beam  a  minute  earlier.  Thus,  from  each  of  the 
16  NORSAR  events,  30  synthetic  events  were  produced  making  a  total  population  of  480  synthetic 
events.  Each  set  of  5  channels  was  bandpass  filtered  with  a  3-pole  Butterworth  filter  with  3-dB 
points  at  0.8  and  3.5  Hz. 

The  event  detector  was  run  on  all  5  channels  for  each  of  the  events.  Figure  VII-4  shows  a 
histogram  of  the  results  as  a  function  of  NORSAR  m^.  The  shaded  area  indicates  the  number  of 
detected  events,  and  the  unshaded  area  denotes  the  undetected  events.  These  data  were  used  to 
compute  the  incremental  detection  probability  for  the  detection  scheme.  First,  the  histogram 
data  were  grouped  in  bins  0.2  magnitude  unit  wide,  then  the  ratio  of  detected  events  to  events 
processed  was  computed  at  each  magnitude.  The  results  are  displayed  in  Fig.  VII-5  where  the 
solid  curve  is  the  detection  probability  for  all  5  channels,  i.e.,  using  the  histogram  of  Fig.  Vll-4, 
and  the  dashed  lines  are  for  the  2  channels  showing  the  largest  deviation  from  the  average  for 
detection  probabilities  >0.5.  For  all  channels  together,  the  90-percent  incremental  detection 
threshold  is  about  =  4,8.  These  results  are  for  a  small  fixed  window  of  winter  noise.  The 
90-percent  threshold  would  be  expected  tc  be  a  little  lower  if  summer  noise  data  were  used. 

An  application  of  this  program  to  the  problem  of  multiple  events  is  being  studied.  In  partic¬ 
ular,  the  constant  updating  of  the  LTA  without  a  lag  after  a  detection  occurs  will  enable  this 
scheme  to  automatically  detect  short-term  bursts  of  energy  in  the  coda  of  other  events. 

C.  W.  Frasier 
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Fig.  VII-1.  Seismicity  of  Central  Asia  in  region  bounded  by  latitudes  18°N  to  58  N 
and  60  °E  to  140  °E,  until  1971.  All  magnitudes  have  mb  >  4.0. 
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Fig.  VII -4.  Histogram  of  events  as  a  function  of  NORSAR  mb.  Shaded  area 
indicates  detected  events,  and  unshaded  area  shows  undetected  events. 


Fig.  VII-5.  Incremental  detection 
probability  for  detection  scheme. 
Data  of  Fig.  VII -4  were  averaged 
over  bins  0.2  magnitude  wide. 
Solid  curve  equals  ratio  of  detec¬ 
tions  to  total  events  using  aver¬ 
aged  data.  Also  shown  are  detec¬ 
tion  probabilities  for  channels 
3  BO  and  4  BO. 
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